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C^ . ABSTRACT 



Based on high-resolution {R « 42 000 to 48 000) and high signal-to-noise (S/N « 50 to 150) spectra obtained with UVESA'LT, we 
present detailed elemental abundances (O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Zn, Y, and Ba) and stellar ages for 12 new microlensed 
dwarf and subgiant stars in the Galactic bulge. Including previous microlensing events, the sample of homogeneously analysed bulge 

i^SC ' dwarfs has now grown to 26. The analysis is based on equivalent width measurements and standard 1-D LTE MARCS model stellar 

atmospheres. We also present NLTE Li abundances based on line synthesis of the 'Li line at 670.8 nm. The results from the 26 
microlensed dwarf and subgiant stars show that the bulge metallicity distribution (MDF) is double-peaked; one peak at [Fe/H] » -0.6 
and one at [Fe/H] a -1-0.3, and with a dearth of stars around solar metallicity. This is in contrast to the MDF derived from red giants 

C ' in Baade's window, which peaks at this exact value. A simple significance test shows that it is extremely unlikely to have such a gap 

^r^ I in the microlensed dwarf star MDF if the dwarf stars are drawn from the giant star MDF. To resolve this issue we discuss several 

possibilities, but we can not settle on a conclusive solution for the observed differences. We further find that the metal-poor bulge 
dwarf stars are predominantly old with ages greater than lOGyr, while the metal-rich bulge dwarf stars show a wide range of ages. 
The metal-poor bulge sample is very similar to the Galactic thick disk in terms of average metallicity, elemental abundance trends, 
and stellar ages. Speculatively, the metal-rich bulge population might be the manifestation of the inner thin disk. If so, the two bulge 

K^ , populations could support the recent findings, based on kinematics, that there are no signatures of a classical bulge and that the Milky 

5_] ■ Way is a pure-disk galaxy. Also, recent claims of a flat IMF in the bulge based on the MDF of giant stars may have to be revised based 

j^ ' on the MDF and abundance trends probed by our microlensed dwarf stars. 
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1. Introduction 
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paper also includes data gathered with the 6^5 m Magellan Telescopes Iq^^^^^Mo^ . Being a major component of nearby galaxies and 

located at the Las Campanas Observatory, Chile, and data obtained at — : T^. ° , ■'. J ^ ^, , .jf , 

»u «7A/r iz 1 <-,u » u- u ■ t A ■ »■« , galaxy populations and a primary feature that classifies galaxies, 

the W.M. Keck Observatory, which is operated as a scientific partner- a j t- t- t- j o > 

ship among the California Institute of Technology, the University of 

California and the National Aeronautics and Space Administration. anonymous ftptocdsarc.u-strasbg.fr (138.79.128.5) or via 

** Table 4 is also available in electronic form at the CDS and http : //cdswe b .u-strasbg . fr/cgi-bin/qcat?J/A-i-A/XXX/AXX| 

full Table 5 is only available in electronic form at the CDS via *** J.A. Johnson is a guest professor at Lund University 



T. Bensby et al.: Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars. IV. 



it is clear that understanding the origin and evolution of bulges 
is integral to the theory of galaxy formation. The central part of 
our own Milky Way harbours a central bulge which enables us 
to study suc h a stellar system in a detail im possible for any other 
galaxy (e.g. jKormendv & Kennicuttl2004l for a review of bulges 
in general). For instance, the next closest bulge stellar system 
is that of the Andromeda galaxy which is more than a hundred 
times more distant. In spite of its "proximity" and the many de- 
tailed spectroscopic and photometric studies during the last few 
decades, the origin and evolutionary history of the Galactic bulge 
is still poorly understood. Its generally very old stellar popula- 
tion, metal-rich nature, and over-abundances of a-elements (e.g., 
Mc William & RichI 119941: I Zoccali et al.1 | 2006t JFulbright et all 






20071: iMelendez et alj|2008[ iBensbv et al.ll2010d) are consistent 
with a classical bulge formed during the collapse of the proto- 
galaxy and subsequent mergers, which w ould have resulted in 
an intense burst of star fo r mation (e.g.. White & Rees 1978^, 
Matteucci & Brocatol 119901: iFerreras et aJ.. .2003.: .Rahimi et al.i 



20101) . Alternatively, the boxy/peanut-like shape of the bulge 



suggests an origin through dynamical instabilities in an already 
establ i shed inner disk ( e.g., Maihara et alj|1978l : ICombes et alj 
Il990t IShen et aTl l2010h . Such secular evolution could possi- 
bly explain the recent discovery of chemical similarities be- 
tween the bulge and t he Galactic thick disk as observed in the 
solar neighbourhood (Melendez et al.l 2008 : Alves-Bri to et al.l 



[2010; Bensby et al. 2010c; Gonzal ez et al J 12011). through the 
action of radial migration of stars dSellwood & BinnevI l2002t 



ISchonrich & Binnevi2009HLoebman et alJl201 ih . 

The metallicity distribution function of the bulge, inferred 
from photometric and spectroscopic studies of red giant stars, 
peaks around the solar value with a significant fraction of super- 
solar metallicity stars and a lo w-metallicity tail extending down 
to at least [Fe/H1«i -1 (e g.. ISadler et all 119961: IZoccali et al] 
l2003Ll2008tlFulbright et al.ll2007h . IndeedTthere are ample indi- 
cations that the Bulge should harbour substantially more metal- 
poor stars (e.g.. lBensbv et al.ii2010b.) . In fact, the very first stars, 
should any have survived to the present day, may well preferen- 
tially be found in t he central regions of Milky Way-type galax- 
ies (e .g., IWvse & G ilmore 1992; Brook etal. 2007; Tumlinson 
I2OIOI) . Unfortunately, the high stellar densities in the bulge and 
its generally high metallicity makes finding such stellar survivors 
from the earliest cosmic epochs very challenging although sev- 
eral large-scale spectroscopic surveys of the Bulge are currently 
unde rway, which may dis cover some of these extremely old stars 
(e.g.. lHowardetalJl2009l) . 

Studies of the detailed chemistry of the bulge have so far 
mainly used intrinsically bright stars (as in the references above). 
However, results based on spectra from giant stars are not trivial 
to interpret as evolutionary processes erase some of the abun- 
dance information. Also, their relatively cool atmospheres result 
in spectra ric h in lines from mole cules, which are difficult to 
analyse (e.g.. lFulbright et ani2006l) . The spectra of dwarf stars, 
on the other hand, even metal-rich ones, are fairly straightfor- 
ward to analy se and are the best trac ers of Galactic chemical 
evolution (e.g. jEdvardsson et al.ll993h . The main difficulty with 
observing dwarf stars in the bulge is their faintness (V = 19 - 20, 
[Feltzing & Gilmore 2000), impeding spectroscopic observations 
under normal circumstances. To obtain a spectrum of high qual- 
ity (R > 40000 and S /N > 75) of a bulge dwarf star would 
require, even if using the largest 8-10 meter telescopes available 
today, an exposure time of more than 50 hours. Microlensing of- 
fers the unique opportunity to observe dwarf stars in the bulge. 
In the event that the bulge dwarf star is lensed by a foreground 
object, its brightness can increase by a factor of several hun- 



dred, making it possible to obtain, in 1 to 2 hours, a spectrum 
of sufficiently high resolution and S /N, adequate for an accurate 
detailed elemental abundance analysis. 

By observing microlensed dwarf stars in the bulge several 
recent studies have giv en a new perspectiv e on t he chemical 
properties of the bulge ( Johnson et al. 200^ 2008h Cohen et al 
2008, 2009, 2010: 'Bensbv et al.l l2009a b. 2010blld:]Epstein et al 
2010) . Main findings so far include the first age-metalUcity 
relation in the bulge which shows that metal-poor stars are 
generally old, but m etal-rich ones have a wide range of ages 
(Bensbv et al. l2010ct) . In addition, abundance ratios for 14 el- 
ements studied in the bulge dwarfs at sub-solar metallicities 
are in excellent agre ement with the abundance patterns in local 
thick disk stars (e.g..|Melendez et al.ll2008l: iBensbv et ani2010cl: 



lAlves-Brito et alJifoToh . Also. iBensbv et all (l2010bl) presented 
the first clear detection of Li in a bul ge dwarf star, showin g that 
the bulge follows the Spite plateau Ispite & Spitdll982l) . The 
most striking result to date from microlensed dwarf stars is that 
the metallicity distribution function (MD F) for dwarf s t ars and 
giant stars in the Galactic bulge diff'ers. In lBensbv et al.l ( 1201 Od) 
we found that the bulge MDF appeared bimodal for the dwarf 
stars, with a paucity at the metallici ty where the MDF b ased on 
giant stars in Baade's window from IZoccali et alJ (l2008h peaks. 
Understanding this discrepancy is vital when studying external 
galaxies where dwarf stars can not be studied, where we have to 
rely on the integrated light from all stars, which is dominated by 
giant stars. 

These results illustrate that observations of dwarf stars pro- 
vide unique information on the evolution of the bulge. For ex- 
ample, the microlensed bulge dwarf stars will have an impor- 
tant impact on the modelling of the bulge, in particular regard- 
ing recent suggestions that the initial mass function (IMF) in 
the bulge needs to be diff'erent from that in the solar neigh- 
bourh ood in order to explain th e MDF based on red giant 
stars JCescutti & Matteuccill201 li) . Additionally, combining the 
dwarf abundances a nd kinematics with numerical studies, e.g., 
iRahimi et al.l (l2010l) . points to the possibility of disentangling 
dilTerent formation scenarios for the bulge, e.g. secular versus 
merger origin. 

In this paper we report the most recent findings from our on- 
going project on the chemical evolution of the Galactic bulge 
as traced by dwarf and subgiant stars that have been observed 
whilst being optically magnified during microlensing events. 
Adding eleven new events, and a re-analysis of MOA-2009- 
BLG-259S, the sample now consists of in total 26 microlensed 
dwarf and subgiant stars in the bulge that have been homoge- 
neously analysed. 



2. Observations and data reduction 

The twelve (11 new and one re-analysis) microlensing events 
presented in this study were all discovered b y the MOA mi- 
crolensing alert svstenrl(e.g.. lBond et al.l200ll) . Most stars were 
observed on the night when the microlensing events peaked, and 
as close to peak brightness (A^ax) as possible (the light curves 
are shown in Fig.[Tli. Exceptions are MOA-2009-BLG-174S and 
MOA-2011-BLG-058S which were observed one day ahead of 
peak brightness as the full Moon would be very close on the 
night of the predicted An,ax- As can be seen in Fig. [Tj the bright- 
nesses for most events were almost constant, or only slowly 
varying, when the spectroscopic observations were carried out. 



https : //itQ19989 . massey . ac . nz/moa/alert/index . html 
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Fig. 1. Light curves for the twelve new microlensing events. Each plot has a zoom window, showing the time intervals when the 
source stars were observed with high-resolution spectrographs. In each plot the un-lensed magnitude of the source star is also given 
(/s). 



The event that showed most variation during the observa- 
tions is MOA-2010-BLG-285S; from / ^ 14.3 at the start of 
the 2 hour exposure to I ^ 13.3 at the end, i.e., it more than 
doubled its brightness during the observations. Furthermore, this 
microlensing event has a binary lens and the UVES observa- 
tions were carried out during , or very close to, a caustic cross- 
ing. In iBensbv et al.l (1201 Obi) , where we presented the first re- 
sults for MOA-201 0-BLG-285S and reported its Li abundance, 
we showed that even though the three UVES exposures not only 
differed in magnification, but also in magnification across the 
stellar surface, the joint effect is an even magnification across 
the stellar surface. Hence, the final co-added spectrum should 
not be affected. As the model light curve for this event requires 
more analysis, which will be published elsewhere, the light curve 
in Fig. [T] currently only contains the observed data points and 
no model fit. Due to a possible planet detection in the lens, the 
model light curve for MOA-2010-BLG-523S also needs more 
analysis, and therefore also only have observed data points in 
Fig.[T] MOA-2011-BLG-034S is an event with a binary lens and 



the microlensing model needs more work. Hence, we have cur- 
rently no A„,ax, Te, or model curve (but only the observed data 
points in Fig.iB for MOA-2011-BLG-2011-034S. 

Ten of the twelve events were observed in 20 10 and 2011 
with the UVES spectrograph (iDekker et al.l l2000l) on the Very 
Large Telescope as part of our ongoing Target of Opportunity 
(ToO) program. They were observed with a 1 .0 arcsec wide slit 
and dichroic #2, giving spectra with a resolution of R ^42 000 
and a spectral coverage between 3760-4980 A (blue ccd), 5680- 
7500 A (lower red ccd), and 7660-9460 A (upper red ccd). A fea- 
tureless spectrum of a rapidly rotating B star (either HR6141 or 
HR7830) was obtained each time. This spectrum was used to 
divide out telluric lines in the bulge star spectrum. The signal- 
to-noise ratios at 6400 A are listed in Tableland range between 
35 and 150. Reductio ns were carried out with the UVES pipeline 
(iBallester et al.l2000l) . version 4.7.8, using the gasgano interface. 

The UVES data were complemented with data for two stars 
observed in 2009. MOA -2009-BLG-259 S was observed with the 
HIRES-R spectrograph (IVogtetal .119941) on the Keck telescopes 
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Table 1. Summary of the twelve microlensing events in the Bulge presented in this study". 



Object 


RAJ2000 


DEJ2000 


/ 


b 


Te 


T 

^ IIIO.\ 


^mcL\ 


Tabs 


Exp. 


SIN 


Spec. 


R 




[hh:mm:ss] 


[dd:mm:ss] 


[deg] 


[deg] 


[days] 


[HJD] 




[MJD] 


[s] 








MOA-2009-BLG-174S 


18:02:07.70 


-31:25:26.40 


-0.34 


-4.34 


65 


4963.8147 


1100 


A^l.-iU 


9000 


95 


MIKE 


42 000 


MOA-2009-BLG-259S 


17:57:57.50 


-29:11:39.00 


1.15 


-2.45 


62 


5016.7674 


204 


5016.395 


1800 


110 


HIRES 


48 000 


MOA-2010-BLG-037S 


18:05:17.93 


-27:56:13.21 


3.04 


-3.23 


92 


5289.4614 


19 


5284.308 


7200 


35 


UVES 


42 000 


MOA-2010-BLG-049S 


18:05:07.02 


-26:46:12.79 


4.04 


-2.63 


12 


5263.7554 


110 


5263.289 


7200 


90 


UVES 


42 000 


MOA-2010-BLG-078S 


17:52:01.65 


-30:24:25.63 


-0.64 


-1.95 


30 


5285.8370 


22 


5285.303 


7200 


40 


UVES 


42 000 


MOA-2010-BLG-167S 


18:11:27.37 


-29:41:02.54 


2.16 


-5.26 


65 


5380.8834 


22 


5380.129 


4740 


100 


UVES 


42 000 


MOA-2010-BLG-285S 


17:56:48.16 


-30:00:39.47 


0.32 


-2.64 


42 


4354.7421 


390 


5354.087 


7200 


150 


UVES 


42 000 


MOA-2010-BLG-311S 


18:08:49.98 


-25:57:04.27 


5.17 


-2.96 


19 


5365.1937 


540 


5365.001 


4800 


55 


UVES 


42 000 


MOA-2010-BLG-446S 


18:07:04.25 


-28:03:57.77 


3.12 


-3.64 


20 


5413.8937 


61 


5413.125 


7200 


70 


UVES 


42 000 


MOA-2010-BLG-523S 


17:57:08.87 


-29:44:58.40 


0.58 


-2.58 


20 


5432.6283 


290 


5431.997 


7200 


150 


UVES 


42 000 


MOA-2011-BLG-034S 


18:10:12.23 


-27:32:02.39 


3.92 


-3.99 


- 


5638.62 


- 


5639.317 


7200 


45 


UVES 


42 000 


MOA-2011-BLG-058S 


18:09:00.36 


-29:14:12.66 


2.30 


-4.57 


17 


5646.92 


90 


5645.320 


7200 


45 


UVES 


42 000 



Notes, t Given for each microlensing event is: RA and DE coordinates (J2000) read from the fits headers of the spectra (the direction where the telescope pointed 
during observation); galactic coordinates (/ and h); duration of the event in days {Te)\ time when maximum magnification occured (Tuiav): maximum magnification 
(A,„fl j); time when the event was observed with high-resolution spectrograph (Tots)', the exposure time (Exp.), the measured signal-to-noise ration per pixel at 
~6400A (SIN); the spectrograph that was used (Spec); the spectral resolution (R). 

Table 2. Stellar parameters and radial velocities for the 12 new microlensed dwarf stars. ' 



Object 


Teff 


logs' 


^. 


[Fe/H] 


A'pel. 


Fell 


V, 




[K] 


[cgs] 


[kms-i] 








[km/s] 


MOA-2010-BLG-285S 


6064 ± 129 


4.20 ± 0.23 


1.85 ±0.38 


-1.23 ±0.09 


53, 


12 


+46.0 


MOA-2010-BLG-078S 


5231 ± 135 


3.60 ±0.31 


1.30 ±0.29 


-0.99 ±0.15 


54, 


3 


+52.3 


MOA-2010-BLG-167S 


5444 ± 58 


4.00 ±0.14 


1.17±0.11 


-0.59 ± 0.05 


118, 


17 


-79.4 


MOA-2010-BLG-049S 


5738 ±58 


4.20 ±0.12 


0.95 ±0.11 


-0.38 ±0.06 


96, 


14 


-116.7 


MOA-2010-BLG-446S 


6327 ± 123 


4.50 ±0.19 


1.57 ±0.17 


-0.37 ± 0.08 


69, 


8 


+56.5 


MOA-2010-BLG-523S 


5250 ± 98 


4.00 ± 0.23 


2.10 ±0.24 


H-0.10±0.07 


54, 


11 


+97.3 


MOA-2009-BLG-174S 


5620 ± 65 


4.40 ±0.12 


1.23 ±0.11 


+0.13 ±0.05 


121, 


17 


-21.4 


MOA-2011-BLG-034S 


5467 ± 67 


4.10 ±0.15 


0.87 ±0.14 


+0.13 ±0.06 


84, 


11 


+ 127.0 


MOA-2009-BLG-259S 


4953 ± 93 


3.40 ± 0.24 


1.02 ±0.11 


+0.37 ± 0.05 


72, 


11 


+81.5 


MOA-2011-BLG-058S 


5260 ± 75 


4.00 ±0.16 


0.67 ±0.14 


+0.39 ± 0.05 


63, 


10 


-139.7 


MOA-2010-BLG-311S 


5464 ± 108 


3.80 ±0.26 


1.29 ±0.19 


+0.49 ±0.11 


58, 


11 


+44.4 


MOA-2010-BLG-037S 


5745 ± 95 


4.00 ± 0.26 


1.51 ±0.20 


+0.56 ±0.10 


56, 


13 


-8.4 



Notes. ' Given for each star is: effective temperature (T^g); surface gravity (log g)\ microturbulence parameter (^t): [Fe/H]; number of Fe i and Fe ii lines used in the 
analysis; radial velicity («r). 



on Hawaii. The slit was 0.86arcsec wide, giving a spectral res- 
olution of /? w 48 000, a wavelength coverage from 3900 to 
8350 A, with small gaps between the orders beyond 6650 A. The 
signal-to-noise ratio is S IN * 110 per pixel aro und 6400 A. 
Detail s of the reduction procedure can be found in ICohen et al.l 
(l2008h . Actu ally, an analys i s of ]V[O A-2009-BLG-259S was first 
presented in iBensbv et al.l (1201 Od) . However, that analysis was 
based on a spectrum that was obtained with UVES and when 
only the UVES blue arm was available due to maintenance work 
on the red arm. That spectrum had a very limited blue wave- 
length coverage and as bulge stars are heavily reddened, it also 
had poor S/N. As a result only a few Fei and Feii could be 
measured and the stellar parame ters were poorly con strained 
with large uncertainties. Hence, in lBensbv et al.l ( l2010d) . IVlOA- 
2009-BLG-259S was excluded when discussing the metallic- 
ity distribution and abundance trends of the microlensed bulge 
dwarf stars. The first met allicity result, base d on this HIRES 
spectrum, was presented in lCohen et al.l(l2010l) . For consistency, 
it has here been re-analysed with the same methods as for all the 
other events. 



MOA-2009-BLG-174S was observed with the MIKE spec- 
trograph (Bernstein et al. 2003) on the Magellan Clay telescope 
at the Las Campanas observatory in Chile. A 0.7 arcsec wide 
slit was used giving a spectrum with R x 4-2 000 and continuous 
coverage between 3500 to 9200 A. The signal-to-noise per pixel 
at 6400 A is around 95. Reductions were made with the Carnegie 
python pipelinqj 

In total we are now in possession of high quality spectra for 
26 microlensed dwarf and subgiant stars in the bulge. Table [1] 
lists the details for the twelve n ew events, while the d ata for the 
other 14 events can be found in lBensbv et al.l(l2010cl) . 



3. Analysis 

3. 1 . Stellar parameters and elemental abundances 

Stellar parameters and elemental abundances have been de- 
termined through exactly the same meth ods as in our previ- 
ous stu dies of microlensed bulge dwarfs ( iBensbv et al.ll2009R 
l2010d) . Briefly, we use standard 1-D plane-parallel model stel- 



Available at |http : //obs . carnegiescience . edu/Code/mike| 
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Table 3. Stellar ages and a comparison of colours and effective temperatures as determined fr om spectroscopy and microlensing 
techniques for all 26 microlensed dwarf stars. The lower part of the table includes the stars from lBensbv et al.l(l2010d) .^ 



Object 



Tefr M logL Age -Icr +lo- (V-I)o My -Icr +lo- M, 
[K] Mo Lq [Gyr] [Gyr] [Gyr] [mag] [mag] [mag] [mag] [mag] 



Mf (V-Ifg T^^ Comments 
[mag] [mag] [K] 



MOA-2010-BLG-285S 


6064 


0.81 


0.04 


11.6 


7.4 


13.3 


0.63 


4.80 


3.78 


5.06 


4.18 


4.84 


0.58 


6283 




MOA-2010-BLG-078S 


5231 


0.84 


0.75 


12.7 


5.0 


15.5 


0.86 


3.04 


2.23 


3.45 


2.18 


1.86 


0.85 


5259 


very high DR 


MOA-2010-BLG-167S 


5444 


0.87 


0.34 


14.0 


8.6 


15.1 


0.79 


4.03 


3.61 


4.11 


3.24 


2.78 


0.72 


5678 




MOA-2010-BLG-049S 


5738 


0.88 


0.16 


11.7 


9.7 


13.3 


0.70 


4.46 


4.10 


4.84 


3.76 


2.74 


0.69 


5785 


some DR 


MOA-2010-BLG-446S 


6327 


1.06 


0.23 


3.2 


1.4 


4.7 


0.56 


4.11 


3.81 


4.45 


3.55 


- 


- 


- 


a bit of DR, corrupt CMD 


MOA-2010-BLG-523S 


5250 


0.99 


0.26 


13.1 


5.5 


13.2 


0.86 


4.29 


3.60 


4.67 


3.43 


- 


0.77 


5510 


a bit of DR 


MOA-2009-BLG-174S 


5620 


0.98 


-0.05 


8.7 


3.8 


9.7 


0.73 


4.99 


4.62 


5.12 


4.26 


3.57 


0.71 


5703 




MOA-2011-BLG-034S 


5467 


0.96 


0.17 


11.6 


8.0 


13.2 


0.77 


4.46 


3.85 


4.86 


3.69 




0.67 


5720 




MOA-2009-BLG-259S 


4953 


1.21 


0.79 


3.0 


1.8 


7.7 


1.00 


3.03 


2.09 


3.89 


2.03 


2.95 


0.82 


5366 




MOA-2011-BLG-058S 


5260 


0.97 


0.34 


12.7 


6.4 


13.1 


0.85 


4.08 


3.77 


4.66 


3.23 


3.36 


0.82 


5366 




MOA-2010-BLG-311S 


5464 


1.09 


0.39 


4.9 


3.1 


9.1 


0.78 


3.92 


2.81 


4.50 


3.14 


3.21 


0.75 


5568 


high DR 


MOA-2010-BLG-037S 


5745 


1.25 


0.49 


4.4 


2.9 


5.5 


0.69 


3.63 


2.98 


4.30 


2.94 


3.61 


0.80 


5421 




OGLE-2009-BLG-076S 


5877 


0.84 


-0.04 


11.6 


7.0 


13.1 


0.67 


5.00 


4.26 


5.17 


4.33 


4.24 


0.70 


5752 




MOA-2009-BLG-493S 


5457 


0.74 


-0.41 


13.1 


4.3 


13.1 


0.79 


5.83 


5.59 


6.11 


5.04 


4.81 


0.70 


5752 


some DR 


MOA-2009-BLG-133S 


5597 


0.81 


-0.24 


12.1 


5.0 


15.9 


0.74 


5.54 


4.70 


5.77 


4.80 


4.24 


0.71 


5715 


double clump sightline 


MOA-2009-BLG-475S 


5843 


0.85 


-0.08 


8.7 


4.5 


12.3 


0.68 


5.10 


4.12 


5.33 


4.42 


4.30 


0.62 


6070 


some DR 


MACHO-1999-BLG-022S 


5650 


0.88 


0.36 


13.1 


7.0 


13.4 


0.73 


3.97 


3.49 


4.70 


3.24 


- 


- 


- 




OGLE-2008-BLG-209S 


5243 


0.93 


0.38 


8.6 


5.3 


12.7 


0.86 


3.96 


3.34 


4.11 


3.10 


2.57 


0.76 


5542 




MOA-2009-BLG-489S 


5643 


0.93 


-0.09 


11.6 


6.1 


12.8 


0.73 


4.97 


4.33 


5.22 


4.24 


3.44 


0.88 


5200 


high DR 


MOA-2009-BLG-456S 


5700 


1.00 


0.05 


8.7 


4.9 


9.9 


0.71 


4.59 


4.00 


4.90 


4.88 


2.81 


0.69 


5773 




OGLE-2007-BLG-514S 


5635 


1.06 


0.16 


6.4 


4.2 


9.6 


0.73 


4.49 


3.53 


4.87 


3.76 


4.41 


0.73 


5634 


high DR, binary peak 


MOA-2008-BLG-311S 


5944 


1.17 


0.18 


2.9 


1.1 


3.5 


0.64 


4.37 


4.06 


4.51 


3.73 


3.98 


0.69 


5767 


a bit of DR 


MOA-2006-BLG-099S 


5741 


1.10 


0.04 


3.3 


1.5 


5.0 


0.70 


4.66 


4.35 


4.87 


3.96 


3.86 


0.77 


5508 


very high DR 


MOA-2008-BLG-310S 


5704 


1.09 


0.07 


4.9 


3.0 


6.3 


0.71 


4.69 


4.04 


4.76 


3.98 


3.51 


0.72 


5664 


Janczak et al. (2010') 


OGLE-2007-BLG-349S 


5229 


0.95 


0.15 


13.1 


9.6 


13.9 


0.87 


4.58 


4.34 


5.30 


3.71 


4.21 


0.81 


5393 


a bit of DR 


OGLE-2006-BLG-265S 


5486 


1.04 


0.10 


8.7 


5.4 


9.5 


0.78 


4.65 


4.23 


5.00 


3.87 


3.64 


0.71 


5696 





Notes. ' Column 2 gives the spectroscopic temperature; col. 3 the stellar mass (inferred from ischrones); col. 4 the luminosity (inferred from isochrones), col. 5 the 
stellar age (inf erred from isochrones); c ols. 6 and 7 the 1-sigma lower and upper age limits; col. 8 the "spectroscopic" colours, based on the colour-[Fe/H]-reff 
calibrations bv lCasagrande et aU I2O10I) : col. 9 the absolute V magnitude based on the spectroscopic stellai' parameters; cols. 10 and 1 1 the 1-sigma lower and upper 
limits on My', col. 12 the absolute / magnitude based on spectroscopic stellar parameters; col. 13 the absolute / magnitude based on microlensing techniques; 
col. 14 the (V — / ) o colo ur based on microlensing techniques; col. 15 the infen'ed effective temperatures using the colour-[Fe/H]-Teff calibrations by 
ICasagrande et all i2O10ll : col. 16 Comments: DR stands for differential reddening. 



Table 5. Measured equivalent widths and calculated elemental 
abundances for the 12 new microlensed dwarfs. ' 



Table 6. _ LTE Li 

iLind et all (120091) . 



abundances and NLTE corrections from 



Element A 

[A] 



[eV] 



star 1 

W,t £(X) 



star 10 

W,l €(X) 



^ For each line we give the log gf value, lower excitation energy (x\), 
equivalent width iW,i), absolute abundance (log e(X)). The table is only 
available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (13Q.79. 128. 5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/XXX/AXX. 



Object 



Iog6(Li)LTE Anlte reff[K] [Fe/H] 



OGLE-2006-BLG-265S' 


0.96 


+0.08 


5486 


+0.47 


OGLE-2007-BLG-349S+ 


0.89 


+0.11 


5229 


+0.42 


]VIOA-2008-BLG-311St 


2.23 


+0.01 


5944 


+0.36 


IVIOA-2010-BLG-049S 


1.43 


+0.07 


5738 


-0.38 


IVIOA-2010-BLG-285S 


2.21 


-0.05 


6064 


-1.23 


IVIOA-2010-BLG-523S 


1.64 


+0.10 


5250 


+0.09 



^ Full set of stellar parameters and elemental abundances for 
OGLE-200 6-BLG-265S, OGLE- 2007-BLG-349S, and MOA-2008-BLG-31 IS 
aie given in lBensbv et all )2010ci) . 



lar atmospheres calculated with the Uppsala MARCS code 
dGustafsson et al.lll975l: lEdvardsson et al.lll993l: lAsplund et alj 
Il997h . (For consistency with our previous anal yses, we continue 
to use the old MARCS models. As shown in iGustafsson et all 
l2008l the differences between the new and old MARCS mod- 
els are very small for our types of stars, i.e. F and G dwarf 
stars.) Elemental abundances are calculated with the Uppsala 
EQWIDTH program using equivalent widths that were mea- 
sured by hand using the IRAF task SPLOT. Effective temper- 
atures were determined from excitation balance of abundances 
from Fei lines, surface gravities from ionisation balance be- 
tween abundances from Fez and Feii lines, and the microtur- 



bulence parameters by requiring that the abundances from Fe i 

lines are independent of line strength. 

An error analysis, as outlined in lEpstein et alJ (1201 Oh . has 
been performed for the microlensed dwarf stars. This method 
takes into account the uncertainties in the four observables that 
were used to find the stellar parameters, i.e. the uncertainty of 
the slope in the graph of Fe i abundances versus lower excita- 
tion potential; the uncertainty of the slope in the graph of Fei 
abundances versus line strength; the uncertainty in the difference 
between Fe i and Fe ii abundances; and the uncertainty in the dif- 
ference between input and output metallicities. The method also 
accounts for abundance spreads (line-to-line scatter) as well as 
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Table 4. Elemental abundance ratios, errors in the abundance ratios, and number of lines used, for the 12 new microlensed dwarf 
stars'. 





[Fe/H] 


[0/Fe]t 


[Na/Fe] 


[Mg/Fe] 


[Al/Fe] 


[Si/Fe] 


[Ca/Fe] 


[Ti/Fe] 


[Cr/Fe] 


[Ni/Fe] 


[Zn/Fe] 


[Y/Fe] 


[Ba/Fe] 


MOA-2010-BLG-285S 


-1.23 


0.55 


-0.06 


0.42 


0.30 


0.30 


0.35 


0.38 


-0.01 


-0.08 


0.28 


0.39 


0.22 


MOA-2010-BLG-078S 


-0.99 


0.65 


-0.02 


0.47 


0.25 


0.35 


0.26 


0.34 


— 


0.00 


— 


— 


0.32 


MOA-2010-BLG-167S 


-0.59 


0.56 


0.07 


0.43 


0.33 


0.22 


0.23 


0.34 


0.07 


0.03 


0.14 


0.07 


-0.18 


MOA-2010-BLG-049S 


-0.38 


0.36 


-0.05 


0.28 


0.28 


0.14 


0.25 


0.31 


0.02 


0.01 


0.23 


0.1 


0.13 


MOA-2010-BLG-446S 


-0.37 


0.21 


0.10 


0.32 


0.23 


0.18 


0.25 


0.34 


— 


0.05 


— 


— 


-0.06 


MOA-2010-BLG-523S 


0.10 


0.25 


0.19 


0.14 


0.33 


-0.01 


0.17 


0.22 


0.16 


0.00 


-0.24 


— 


-0.18 


MOA-2009-BLG-174S 


0.13 


-0.03 


-0.03 


0.09 


-0.04 


0.03 


0.02 


0.08 


0.04 


-0.03 


-0.05 


-0.03 


-0.04 


MOA-2011-BLG-034S 


0.13 


-0. 1 1 


0.02 


0.03 


0.07 


0.07 


-0.01 


0.05 


-0.02 


0.03 


0.03 


— 


-0.14 


MOA-2009-BLG-259S 


0.37 


0.03 


0.42 


0.36 


0.35 


0.10 


0.17 


0.14 


0.11 


0.18 


0.41 


0.38 


-0.07 


MOA-2011-BLG-058S 


0.39 


-0.12 


0.20 


0.11 


0.05 


0.08 


0.00 


0.00 


0.10 


0.13 


0.28 


- 


-0.12 


MOA-2010-BLG-311S 


0.49 


-0.24 


0.26 


0.24 


0.18 


0.08 


0.13 


0.07 


0.03 


0.14 


0.14 


— 


-0.13 


MOA-2010-BLG-037S 


0.56 


-0.40 


0.31 


0.16 


0.10 


0.08 


-0.04 


0.14 


0.06 


0.13 


-0.23 


— 


-0.10 




0-[Fe/H] 


0"[0/Fe] 


0"[Na/Fc] 


0"[Mg/Fcl 


0"[A1/Fcl 


0-[Si/FcJ 


O-lCa/Fc] 


O" [Ti/Fe] 


a"[Cr/Fe] 


0"[Ni/FcJ 


a"[Zn/FcJ 


0-[Y/Fe] 


''■[Ba/Fe] 


MOA-2010-BLG-285S 


0.095 


0.26 


0.05 


0.07 


0.07 


0.07 


0.06 


0.12 


0.19 


0.11 


0.10 


0.24 


0.21 


MOA-2010-BLG-078S 


0.15 


0.49 


0.10 


0.13 


0.12 


0.17 


0.19 


0.05 


— 


0.15 


— 


— 


0.14 


MOA-2010-BLG-167S 


0.05 


0.19 


0.09 


0.08 


0.04 


0.06 


0.09 


0.07 


0.07 


0.06 


0.14 


0.20 


0.09 


MOA-2010-BLG-049S 


0.06 


0.18 


0.15 


0.06 


0.03 


0.05 


0.10 


0.09 


0.09 


0.06 


0.10 


0.21 


0.08 


MOA-2010-BLG-446S 


0.08 


0.18 


0.05 


0.06 


0.07 


0.05 


0.07 


0.03 


— 


0.09 


— 


— 


0.24 


MOA-2010-BLG-523S 


0.07 


0.22 


0.13 


0.12 


0.10 


0.08 


0.15 


0.11 


0.08 


0.07 


0.13 


— 


0.11 


MOA-2009-BLG-174S 


0.05 


0.15 


0.08 


0.07 


0.05 


0.04 


0.09 


0.07 


0.05 


0.05 


0.12 


0.26 


0.11 


MOA-2011-BLG-034S 


0.06 


0.19 


0.03 


0.09 


0.06 


0.06 


0.10 


0.07 


0.04 


0.06 


0.14 


— 


0.10 


MOA-2009-BLG-259S 


0.05 


0.24 


0.22 


0.19 


0.10 


0.08 


0.20 


0.09 


0.07 


0.06 


0.14 


0.57 


0.20 


MOA-2011-BLG-058S 


0.04 


0.18 


0.09 


0.09 


0.10 


0.06 


0.10 


0.10 


0.08 


0.05 


0.09 


- 


0.05 


MOA-2010-BLG-311S 


0.11 


0.27 


0.12 


0.19 


0.14 


0.09 


0.16 


0.09 


0.08 


0.09 


0.17 


— 


0.28 


MOA-2010-BLG-037S 


0.10 


0.29 


0.16 


0.18 


0.11 


0.07 


0.12 


0.10 


0.12 


0.09 


0.24 


— 


0.28 




A'Fe, 


No 


^Na 


Nug 


Nm 


A's, 


Wca 


Nj, 


A'Cr 


Nm 


Nzn 


A'y 


Nb. 


MOA-2010-BLG-285S 


53 


3 


1 


6 


3 


13 


10 


15 


3 


11 


2 


3 


4 


MOA-2()10-BLG-()78S 


54 


2 


1 


4 


4 


16 


10 


2 


— 


9 


— 


— 


2 


MOA-2010-BLG-167S 


118 


3 


2 


6 


6 


23 


15 


25 


9 


33 


3 


3 


4 


MOA-2010-BLG-049S 


96 


3 


2 


7 


6 


25 


13 


13 


2 


34 


3 


2 


3 


MOA-2010-BLG-446S 


69 


3 


2 


5 


5 


19 


10 


3 


— 


15 


— 


— 


3 


MOA-2010-BLG-523S 


54 


3 


2 


4 


6 


20 


11 


2 


5 


25 




— 


3 


MOA-2009-BLG-174S 


122 


3 


4 


4 


7 


27 


18 


16 


12 


34 




3 


3 


MOA-2011-BLG-034S 


84 


3 


2 


5 


6 


23 


11 


6 


2 


30 




— 


3 


MOA-2009-BLG-259S 


64 


3 


3 


3 


3 


17 


11 


13 


13 


29 




2 


2 


MOA-2011-BLG-058S 


63 


2 


2 


5 


6 


24 


12 


6 


5 


30 




- 


3 


MOA-2010-BLG-311S 


58 


3 


2 


5 


6 


27 


10 


6 


5 


26 




— 


3 


MOA-2()10-BLG-()37S 


56 


3 


2 


5 


6 


25 


11 


6 


3 


25 




— 


3 




itios for oxygen have been corrected for NLTE effects according to the empirical formula givf 










Note that the abundance r 


;n inlBensbv et all (2004bl). 





how the average abundances for each element reacts to ch anges 
in the stellar parameters. Compared to lEpstein et al.l ( 1201 Ol) . who 
used Ca i lines to constrain the microturbulence parameter, we 
use the same Fe i lines that were used for the effective tempera- 
ture. Although the variables are not fully independent, they are 
only weakly correlated. Thus the total error bar is dominated by 
the uncertainties in the slopes. 

Stellar parameters and error estimates are given in Table |2] 
elemental abundance ratios in Table |4] and measured equivalent 
widths and abundances for individual lines in Table |5] 

We have also determined NLTE Li abundances in 6 of the 26 
microlensed dwarf stars through line profile fitting of the ^Li i 
resonance dou blet line at 670.8 n m. The 1-D NLTE corrections 
are taken from iLind et al.l (l2009l) . The Li line was not detected 
in the remaining 20 stars, and the spectra are not of sufficient 
quality to estimate interesting upper limits. The metho d we use 
to deter mine Li abundances is fully described in Bensbv et al.l 
(1201 Obi) where we presented the Li abundance in MOA-2010- 
BLG-285S. The Li abundances are listed in Table|6] 



3.2. Stellar ages and absolute magnitudes 

The determination of stellar ages has been greatly improved in 
relation to our previous work, by the use of probability dis- 
tribution functions. Fu rthermore, the Yonsei-Yale isochrones 
dPemarque et al.l 120041) have been calibrated to reproduce the 
solar age and solar mass for the input solar stellar parameters 
(Feff/log g/[Fe/H] = 5777 /i:/4.44/0.00). This is described in 
detail in a forthcoming paper (Melendez et al., in prep.). In short, 
for a given set of "observed" stellar parameters (Teff, log g, 
[Fe/H]) and theoretical isochrone points (effective temperature 
T, logarithm of surface gravity G, and metallicity M), the age 
is determined from an isochrone age probability distribution 
(APD): 



dP(age) 



1 



A(age) 



2 p(T,s,\ogg,[Fe/n],T,G,M). 



(1) 



A(age) 



where A(age) is an adopted step in age from the grid of 
isochrones and: 



p <x exp[-(reff - Ty-/2{AT,fff] x 
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Fig. 2. The 12 new microlensed stars plotted on the a-enhanced isochrones from iDemarque et alj ( |2004|) . Each set of isochrones 
have been calculated with the same metallicity and a-enhancement as derived for the stars. In each plot the solid lines represent 
isochrones with ages of 5, 10, and 15 Gyr (from left to right). Dotted lines are isochrones in steps of 1 Gyr, ranging from 0.1 Gyr 
to 20 Gyr. EiTorbars represent the uncertainties in T^s and loggf as given in Table |2l Note that ages reported in Table [3] need not 
coincide exactly with the ages that may be read off directly from the figures. The age determinations are based on probability 
distribution functions and is described in Sect. [3] 



exp[-{logg -Gf/2(Alog gr]x 
exp[-([Fe/H] - M)V2(A[Fe/H])^] . 



(2) 



The errors in the derived stellar parameters (ATf^g, etc.) are 
listed in Table |2l The sum in Eq. ([T) is made over a range of 
isochrone ages and in principle all values of T, G, M. In prac- 
tice, the contribution to the sum from isochrone points farther 
away than Jeff ± 3Areff, etc., is negligible. Therefore, the sum 
is limited to isochrone points within a radius of three times the 
errors around the observed stellar parameters. The probability 
distributions are normalised so that Y^dP - 1 . The most prob- 
able age is obtained from the peak of the distribution while 1 cr 
limits are derived from the shape of the probability distribution. 
Similar formalisms allow us to infer the stellar mass, absolute 
magnitudes, and luminosities. 

For three very old stars (> 12 Gyr), both the peak and the 
upper error bar of the APD are not well defined. In those cases 



we adopted as the age the average of the APD peak and the me- 
dian of the solutions, and the upper error bar was adopted as the 
scatter of the solutions. The probability distributions for mass, 
absolute magnitudes, and luminosities are all well behaved. 

Using the new improved age determination method we have 
re-determined the ages for all 26 stars. The ages, absolute magni- 
tudes, and luminosities are reported in Table[3j and Fig.|2]shows 
the twelve new microlensed dwarf stars together with the rele- 
vant Y2 isochrones. As can be seen, a majority of the stars are 
either main sequence turn-off" stars or subgiant stars, while two 
or three are on the verge of ascending the giant branch. The latter 
ones have just left the subgiant branch, and should be un-affected 
by the internal mixi ng processes that occ ur further up on the gi- 
ant branch (see, e.g.. iPinsonneaultJl 1 997l for a review on mixing 
in stars). 

The masses and luminosities presented in Table [3] will yield 
a "physical" log g, computed from the fundamental relationship 
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Fig. 3. Synthesis of the Ho- line at 6563 A for the three stars that show large deviations between the spectroscopic effective tem- 
peratures and the effective temperatures based on microlensing techniques. Observed spectra are shown with black lines, and the 
synthetic spectra with blue Unes. Left-hand side shows the synthetic profiles based on spectroscopic temperatures and the right-hand 
side the profiles based on temperatures from microlensing techniques. The microlensing temperature for MOA-2010-BLG-311S is 
the one based on the original {V-I\) colour (i.e. T^ff = 4750 K). Note that no attempts to match individual lines has been made. The 
compo sition of the model atmospheres are simply scaled with [Fe/H] relative to the standard solar composition by iGrevesse et al.l 
(l20a7h . 
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Fig. 4. Comparison between out spectroscopic surface gravities 
and the "physical" loggf that can be computed from the lumi- 
nosities and stellar masses inferred from isochrones. On aver- 
age the physical log g;s are 0.03 dex higher (with a dispersion of 
0.09 dex). 



that relates lumino sity, stellar mass, a nd eff'ective temperature 
(see, e.g., Eq. (4) in lBensbv et al.ll2003h . Because a star with the 
best-fit age may not Ue exactly at the spectroscopically measured 
Tetf and log g, the the physical and spectroscopic log g:s agree 
within 0.03 + 0.09 dex, rather than exactly, see Fig.|4] 

Probabilistic age determinations like ours are known to 
suffer from systematic biases mainly related to the sampling 



of isochrone data points (e.g., iNordstrom et al.l |2004 their 
Sect. 4.5.4). Bayesian methods have been developed to tackle 
these problems (e.g..|Pont & Ever 2004 ; J0rgensen & LindegrenI 
I2OO5I: ICasagrande et al.l l201lh but we do not follow this ap- 
proach and therefore our ages could be in principle further im- 
proved. Nevertheless, our current age determination is consistent 
with our previous work on thin/thick disk stars so, in a relative 
sense, our discussions regarding the ages of our sample stars are 
still reliable. 



3.3. Effective temperatures from microlensing techniques 



As in iBensbv et al.l (1201 Od) we confront the effective tempera- 
tures derived from spectrosco pic princ iples with those deduced 
from microlensing techniques dYoo et al. 2004). The microlens- 
ing method assumes that the reddening towards the microlensed 
source is the same as the reddening towards the r ed clump, and 
that the red clump in the Bulge has (V-I)o - 1 .08 (iBensbv et al.l 
l2010cl) and M/ = -0.15 (David Nataf, private communica- 
tion). The de-reddened magnitude and colour are then derived 
from the off'sets between the microlensing source and the red 
clump, in the instrumental colour-magnitude diagram (CMD). 
As the lenses for the MOA-2010-BLG-285S and -523S events 
likely are binaries, and need more work, no absolute magni- 
tudes could be estimated for these stars. Also, the {V - /)o 
colour for MOA-2010-BLG-446S could not be determined be- 
cause there are too few clump stars close to this star in the 
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Fig. 5. (a) Difference between the spectroscopic temperatures 
and the temperatures from microlensing techniques versus ef- 
fective temperature, and (b) the difference between {V - I)o 
colours determined from microlensing techniques and those that 
can be determined from the spectroscopic temp eratures using 
the IRFM calibration by iCasagr ande et al] (l201d) . The average 
differences are {AT^ff} = -80K with a dispersion of 217K, 
and <A(y - /)o> = 0.026 with a dispersion of 0.073 mag. The 
three stars that show largest deviations between microlensing 
and spectroscopic parameters are marked. For MOA-2010-BLG- 
3 11 S (mb 1 03 11 ) the original values are shown with open circles, 
and the corrected value with filled circles (connected with the 
original values by a dash-dotted line). 



colour- magnitude diagram . From the colour- [Fe/Hl-Teff calibra- 
tions by ICasagrande et al] ( 1201 Ob we can determine which tem- 
peratures the (V - /)o colours correspond to. The absolute de- 
reddened magnitudes and colours and the effective temperatures 
are given in Table [3] 

Even though the effective temperatures from microlensing 
techniques and spectroscopic principles are in reasonable agree- 
ment, there are three events for which the differences are dis- 
turbingly large; more than 400 K for MOA-2009-BLG-259S and 
MOA-2009-BLG-489S, and more than TOOK for MOA-2010- 
BLG-311S (Fig. |5ji). Note that stel lar parameter s for M OA- 
2009-BLG-489S were published in iBensbv et al] ( l2010cl) . To 
check the temperatures for these stars, and get a third opin- 
ion, we calculate synthetic spectra for the Ha line wing pro- 
files with the latest versio n (v2.1 oct 2010) of the Sp ectroscopy 
Made Easy (SME) tool (IValenti & Piskunovl \l996ii using the 
latest 1-D LTE plane-parallel MARCS model stellar atmo- 
spheres (Gustafsson et al. 2008i), and line data from VALD 
dPiskunov et al.lll995HKupka et allll999h . In Fig. [3] we show the 
observed line profiles and synthetic line profiles for the Balmer 
Ha line for these three stars, using both the spectroscopic tem- 
peratures (left-hand side plots) and the temperatures from mi- 
crolensing techniques (right-hand side plots). We see that the 
synthetic spectra based on spectroscopic temperatures reproduce 
the wings of the Ha line almost perfectly, while the spectra based 



on temperatures from microlensing techniques clearly do not 
match the observed wing profiles. 

We could not find anything that was obviously wrong with 
the photometry of MOA-2009-BLG259S and MOA-2009-BLG- 
489S, although differential reddening is noticeable in the field of 
MOA-2009-BLG-489S (Nataf, private communication, see also 
Table |3]l. However, for MOA-2010-BLG-311S there was only 
one highly-magnified V point, so the {V-I)o colour co uld not be 
robust ly estimated by the standard techniques used bv lYoo et al.l 
(l2004h . In an alternative attempt to determine the (V -I)o colour 
for MOA-2010-BLG-31 IS we first determined the instrumental 
(I-H) colour, making use of the //-band data that are simultane- 
ously taken with / band on the SMARTS ANDICAM camera at 
the 1.3 m CTIO telescope. We then converted from instrumental 
(/-//) to (V-I) by constructing a VIH colour-colour diagram for 
field stars, which is analogous to the technique for optical colour 
transformations (iGould et al.ll2010l) . Now we find a colour of 
{V - /)o = 0.75 for MOA-2010-BLG-311S, which is in good 
agreement with the "spectroscopic" value of 0.78. Note, how- 
ever, that this method cannot be blindly applied to VIH transfor- 
mations because (at the depth of the SMARTS //-band images) 
the colour-colour diagram is constructed from giants and sub- 
giants, while the transformation applies to dwarfs. Additionally, 
there is a giant-dwarf bifu rcation of the VIH co lour-colour re- 
lations for (y - K)o > 3 (lBessell&B rett"l988'). which corre- 
sponds to (V - /)o > 1.3 (see, in particular, Dong et al. 200(3). 
Nevertheless, since most microlensed dwarfs (including this 
one), have bluer colours than this limit, the method can usually 
be applied, if necessary. 

For the 24 microlensed dwarf stars that have microlens- 
ing (V - /)o colours, including the revised colour for MOA- 
2010-BLG31 IS, we find that our spectroscopic temperatures are 
80 + 217 K lower than the ones based on the colour- [Fe/H]- 
Teff relationships (see Fig. [5^). Note that th e colou rs for the 
previous microlensing events in iBensbv et al.l (1201 Od) were es- 
timated under the assumption that the red clump in the bulge 
has (V - /)o = 1.05. Hence, those colours have been adjusted 
(by adding 0.03 mag) and new effective temperatures were de- 
termined. For cla r ity, al l 26 events are listed in Table |3] The 
ICasagrande et al.l (1201 Oh calibrations can also be used to see 
what (V - /)o colours the spectroscopic effective temperatures 
correspond to. For the new stars these are listed in Table [3] 
and the comparison between photometric and "spectroscopic" 
{V - /)() colours are shown in Fig. |5] On average the spectro- 
scopic colours are in good agreement with the ones from mi- 
crolensing techniques being only slightly lower (the difference 
is 0.02 + 0.07 mag, see Fig.|5]3). 

This d ifference is due t o that we no w use the calibrations 
by ICasagrande et a ll (|201 0|), while we in IBensbv et alj (1201 Pel) 
used th e ones by [R amirez & Melendej (|2005|) . If we were to 
use the iRamirez & Melendea (l2005l) calibrations the difference 
would be 0.00 + 0.06 mag. This means that the {V - /)o value 
for the bulge red clu mp that was revised from 1.05 to 1.08 in 
IBensbv et al.l (1201 Od) could be revised back again to 1 .06. 



4. Results 

4. 1 . The metallicity distribution 

In Fig. |6^ we show the MDF for our sample of 26 dwarf stars 
in the bulge. Since the typical error in [Fe/H] is about 0.1 dex, 
we choose a bin size for the histogram of 0.2 dex. This ensures 
that any feature that we detect in the distribution is not due to 
error statistics, i.e., features that could originate from the mea- 
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Fig. 6. a) MDF for our 26 microlensed dwarf star s in the bulge. 
b) MD F for 204 RGB stars in the bulge, from IZoccali et all 
(l2008h . c) The significance for the failure of reproduce a gap 
vs. position of a gap in the dwarf star MDF (see explanation in 
Sect. mil. 
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Fig. 7. Top panel show ages versus [Fe/H] for the now in to- 
tal 26 microlensed dwarf and subgiant stars in the Bulge. Stars 
with [Fe/H] > are marked with black circles and stars with 
[Fe/H] < with grey circles. Bottom panel shows age his- 
tograms for the metal-rich stars (hatched histogram) and metal- 
poor stars (grey histogram). 



surement uncertainties are suppressed by the size of each bin. We 
note that this MDF appears bimodal with peaks at [Fe/H] a; -0.6 
and K +0.3 dex, and a gap of only two stars at [Fe/H] x; -0.1. In 
Fig.|6j5, for co mparison, we show th e MDF for 204 RGB stars in 
the bulge from lZoccah et al.l (|2008|) . Given that the MDF for the 
RGB stars only has one prominent peak (although asymmetric), 
there is a significant difiference between the RGB and the dwarf 
star MDFs. A Kolmogorov-Smirnov test between the two dis- 
tributions yields a p-value of 0. 11 . Even though this is not low 
enough to reject the null-hypothesis that the MDFs are drawn 
from the same distribution, under the commonly used limit of 
0.05, it is low enough to warrant further investigation. 

The most prominent difference is the deficiency of stars in 
the dwarf star MDF, situated almost where the RGB MDF peaks. 
The question we ask is; what is the probability of obtaining such 
a deficiency, similar to the one in the dwarf star MDF, if 26 stars 
are drawn randomly from the RGB MDF? To answer this ques- 
tion, we construct a simple significance test. First, we draw 26 
stars randomly from the RGB MDF. Given that there are 204 
stars in the RGB sample, we assume the RGB MDF to be rep- 
resentative of the complete underlying MDF. Next, we check if 
the sample of 26 randomly drawn stars has a deficiency of stars 
in the MDF at [Fe/H] = x. We define the deficiency of stars in 
the MDF as a gap with two stars or less in a metallicity range 
of 0.3 dex (see Fig. 5a). The gap in Fig. |6] spans about 0.4 dex. 
However, given that the typical eiTor in [Fe/H] for the dwarf stars 
is about 0.1 dex, we use a range of 0.3 dex. We iterate this pro- 
cess 10^ times. We define the significance as the fraction of it- 
erations that fail to reproduce a gap in the randomly generated 
dwarf star MDF. Finally, we sample over the position of the gap, 
X, from [Fe/H] - -2 to [Fe/H] = +1 in steps of 0.1 dex. 

Figure 5c shows our result from the test. As can be seen, the 
existence of a gap of stars in the dwarf star MDF at [Fe/H] = 



-0.1 is unlikely. This indicates that the gap in the MDF is most 
likely not caused by low number statistics. Additionally, we ran 
the test with a gap range of 0.2 and 0.4 dex and found that it 
did not afi'ect the result significantly. Given the possibility of 
systematic offsets in [Fe/H] for giants and dwarfs, how would a 
shift in the giant MDF affect the significance level as a function 
of [Fe/H]? As we sample over all [Fe/H] in the test above, the 
only effect of shifting the giant MDF is a shift of the significance 
curve in Fig.|6]by the same amount. 

4.2. Ages and metallicities 

Figure |7^ shows the age-metallicity diagram for the 26 mi- 
crolensed dwarf and subgiant stars. At sub-solar metallicities 
the stars are predominantly old with ages between 9 and 13 Gyr 
(see Fig. [Tjj), and no apparent trend of increasing or decreasing 
ages with [Fe/H] is seen. The average age and Icr dispersion is 
11.0 + 2.9 Gyr. This is similar to what is found for thick disk 
dwarfs at the same metallicity in the solar neighbourhood (e.g., 
lFuhrmann|[T99llBensbv et al.ll2007h . 

Furthermore, looking again at the metal-poor population, all 
stars, except MOA-2010-BLG-446S at [Fe/H] = -0.37, seem to 
have the same age within the eiTorbars. If we exclude this star 
the average age becomes 1 1 .7 + 1 .7 Gyr, indicating that the age 
spread in metal-poor population is indeed very small. Regarding 
MOA-2010-BLG-446S, we do not rule out the possibifity that it 
is an interloper from the Galactic disk. 

The stars at super-solar metallicities on the other hand show 
a wide range of ages from the youngest disk, being only a few 
billion years old, to the oldest halo, as old as the Universe (see 
Fig. 17}. The average age and Icr dispersion is 7.6 + 3.9 Gyr for 
the stars at super-solar [Fe/H]. 
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If the ages are weighted by the inverse squares of their errors, 
the average ages and 1 cr dispersions become 5 .9 + 3 .4 Gyr for the 
metal-rich stars, and 11.7 ± 1.5 for the metal-poor stars (9.5 + 
4.2 Gyr when including MOA-2010-BLG-446S). 

4.3. Absolute magnitudes 

In order to further investigate whether our microlensed dwarf 
and subgiant stars are indeed in the Bulge, we have compared 
the absolute / magnitudes obtained from the microlensing pho- 
tometry with the absolute magnitudes obtained from the spectro- 
scopic stellar parameters. 

The spectroscopic absolute magnitudes were obtained with 
the same methods and the same set of isochrones employed to 
determine the ages (see Sect. O. The absolute magnitudes from 
the microlensing technique, M/,^fe„j, were obtained by compar- 
ing the instrumental / magnitude of the microlensing source, 
Iinstr,s, with the instrumental / magnitude of the bulge red clump, 
Iinstr,ciump, assuming that the absolute magnitude of the clump is 
Miciump = -0.15, i.e. 



M, 



l,filens 



= ii, 



insTi;s ' in sfr, clump 



-0.15 



(3) 



This means that the Mj^^ums will exactly equal M^ 



,fi source 



if: 



1. The instrumental clump centroid is correctly measured. 
[Clump centroids are probably accurate to about 0.1 mag. 
This is the typical measured difference in the clump cen- 
troids of different data sets of the same field and then trans- 
formed from one to another ] 

2. The instrumental source magnitude is correctly measured. 
[In most cases this is extremely accurate (of order 0.02 mag 
or less).] 

3. The extinction towards the source is the same as the mean 
extinction towards the clump. [Most fields have modest dif- 
ferential extinction. There is very little systematic bias. The 
random error due to differential extinction is probably less 
than 0.05 mag.] 

4. The absolute magnitude of the bulge clump is Mj^ciump - 
-0.15. [Note; Mi^ciump - -0.25 is the value that usually has 
been used in bulge microlensing studies. However, recent 
studies indicate that it probably is incorrect, and should be 
-0.15 (D. Nataf, private communication).] 

5. The source is at the same distance as the clump. [The sources 
are certainly not at the same distance as the clump. They 
are thought to be in the bulge, and therefore at a range of 
distances similar to the range of distances of clump stars. 
However, the sources must be behind the lenses. Thus, if 
the latter are in the bulge (thought to be about 60% of lens- 
ing events) then the former must be preferentially towards 
the back of the bulge. Hence, on these grounds alone one 
would expect them to be of order 0. 1 mag behind the clump. 
However, from Fig.[TT]it is clear that the spectroscopic sam- 
ple is drawn preferentially from positive longitudes, which 
constitute the near side of the Galactic bar, which has a slope 
of roughly 0.1 mag deg"'. We take this as indicating a strong 
bias towards brighter sources (which then require less mag- 
nification to enter our sample). Of course, within any given 
line of sight, this bias is working on a much narrower base- 
line, i.e., the intrinsic dispersion in the bulge depth, which 
is order 0. 15 mag. We therefore estimate that the two biases 
approximately cancel out, and that there is an intrinsic dis- 
persion of source distances of 0. 1 5 mag. ] 

Figure|8]shows the difference between the absolute / magni- 
tudes from the two methods versus stellar mass (as determined 
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Fig. 8. The difference between the absolute / magnitudes from 
microlensing techniques and from spectroscopy versus stellar 
mass (derived from spectroscopy). Error bars represent the un- 
certainties in the spectroscopic values. 



from spectroscopy). The spectroscopic values are on average 
higher by -0.13 mag (with a dispersion of 0.56 mag). The error 
bars represent the uncertainties in the spectroscopic values. Even 
though the individual points appears to be basically consistent at 
the 1 -cr level, there is a tentative trend in the sense that the spec- 
troscopic values are higher for low masses while microlensing 
values are higher for high masses. 

To formally evaluate the significance of the difference be- 
tween the microlensing and spectroscopic estimates of M/, we 
first estimate an error of 0.2 mag in the microlensing M/ values, 
which is the quadrature sum of the 5 enumerated error sources 
identified above. We add this in quadrature to the (asymmetric) 
errors shown in Fig. [8] and evaluate ;^f^ as a function of assumed 
offset. We find AM, = -0.10 ± 0.15, with y^ . = 12.48 for 21 

* ' ^ mm 

degrees of freedom. The fact that the mean is within 1 cr of zero 
implies that the sample is consistent with being entirely drawn 
from bulge stars. Since x^ is less than the number of dof, there 
is no evidence for scatter beyond that due to measurement er- 
rors (which includes our estimate of a 0.15 mag dispersion in 
distances relative to the clump). Of course, this does not prove 
that all of our sources are in bulge; there could be disk inter- 
lopers. However, the fact that microlensing models place the 
overwhelming majority of lensed sources in the bulge, combined 
with the demonstrated consistency between the M/ determined 
from spectra and those determined from microlensing (based on 
the assumption that they have the same mean distance as the 
bulge clump) strongly supports a bulge location for the great ma- 
jority of our sources. 

4.4. Colour-magnitude diagram 

Figure |9] shows our stars in the HR diagram. We consider the 
metal-poor and metal-rich stars separately, and we include de- 
reddened colours and absolute magnitudes estimated directly 
from the spectra (filled circles, see Sect. 13.2b as well as from 
the micro-lensing technique (open circles, see Sect. I3.3l l. The 
two values are connected for each star by a dotted line (note 
that 4 stars do not have microlensing values, see Table [3]l. The 
metal-poor stars (Fig |9^) are spread from the main-sequence to 
the bottom of the red giant branch. The age of the population of 
these stars is younger if we consider the magnitudes and colours 
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Fig. 9. (a) and (b) show the {V - I)o versus M/ colour-magnitude 
diagram using values from microlensing techniques (open cir- 
cles) and spectroscopy (filled circles), (a) shows the stars 
with [Fe/H] < and (b) the stars with [Fe/H] > 0. Each 
CMP also show the Y- is ochrones (1, 5, 10, and 15Gyr) 
from iDemarque et alJ (l2004l) . For the metal-poor sample the 
isochrones have [Fe/H] - -0.6 and for the metal-rich sample 
[Fe/H] = +0.3. The spectroscopic and microlensing values have 
been connected with a dashed fine for each star (4 stars do not 
have microlensing values). 



derived by micro-lensing techniques (between 5 and 10 Gyr), 
while, if we use the spectroscopically determined colours and 
magnitudes, the age is about 5 Gyr larger with a very nicely de- 
fined sub-giant branch. Regradless of these differences it is clear 
that we are dealing with a population that is old. 

For the metal-rich stars (Fig |9j5) the situation is somewhat 
dififerent. Here we have no stars that have evolved sufficiently 
to trace the sub-giant branch. Additionally, the spectroscopy and 
micro-lensing techniques yield rather similar results for almost 
all stars. The stars are clustered in the turn-off region and there 
seems to be both rather young stars (5 Gyr with very little room 
for a large change) up to stars as old as 10 Gyr. There are quite 
a few stars that potentially cluster around the 1 Gyr isochrone as 
well. 

In summary, when we consider the high and low metallicity 
stars as two stellar populations, and study their location in an HR 
diagram, we see that the metal-poor population is fairly evolved 
and includes essentially older stars, whilst the metal-rich popu- 
lation is less evolved and tentatively includes some rather young 
stars. 



4.5. Abundance trends 

4.5.1 . Similarities to tine Galactic thick disk 

From the analysis of OGLE-2008-BLG-209S, which was the 
first microlensed subgiant star with a metallicity below solar, 
we got the first hint from a chemically un-evolved star that 
the metal-poor parts of the bulge might be s imil ar to those of 
the Gal actic thick disk (iBensbv et al.ll2009bh . In iBensbv et all 
(1201 Od) where the sample was expanded to 15 microlensed 
dwarf and subgiant stars, out of which 7 had [Fe/H] < 0, it 
was clear that the metal-poor bulge, in terms of both elemental 
abundance ratios as well as stellar ages, is similar to the Galactic 
thick disk. With this study we add another five stars with sub- 
solar metallicities, and another seven with super-solar metallic- 
ities. The updated abundance plots are shown in Fig. [TO] and it 
is clear that the new stars further strengthens the similarity be- 
tween the metal-poor bulge and the Galactic thick disk as seen in 



the solar neighbourhood. The agreement is almost perfect for all 
the a-elements and only for the two .s-process elements Ba and 
Y are there some discrepancies, and then for only one or two 
stars at the lowest [Fe/H]. Given the error bars and the fact that 
it is only 1-2 stars, it is difficult to judge if the differences for 
Ba and Y are real. We also note the large dispersion in [Zn/Fe] 
at high metallicities, which most likely is caused by the fact that 
the Zn abundances for a majority of the stars are based on only 
one spectral line. 

The similarity between the bulge and the local thick disk 
abundanc e trends was actually first recognise d from studies of 
K gian ts ('Melendez et al.' 2008) and later by lAlves-Brito et all 
(1201 Ol) and Gonzalez et al. (2011) . Recently it has been shown 
that this similarity also holds when comparing bulge giant stars 
with K giant stars fro m the thick disk in the inner parts of the 
Galaxy (IBensbv et a02010al) . 



4.5.2. Rapid enrichment and time scale for SNIa in the bulge 

The high level of a-element abundances relative to Fe are signa- 
tures of rapid star formation where massive stars, which have 
short lifetimes, are the main contributors of chemical enrich- 
ment of the interstellar medium when they explode as core col- 
laps e SNe produce and ex pel large amounts of cf-elements (see, 
e.g.. lBallero et alj|2007bl) . This will result in high [a/Fe] ratios 
at low [Fe/H], and we see this for the microlensed dwarf stars 
at [Fe/H] < -0.4. Once low- and intermediate-mass stars starts 
to play a significant role in the chemical enrichment, through 
the explosions of SNIa, which produce relatively little of the a- 
elements, there will be a downturn (knee) in the [a'/Fe]-[Fe/H] 
abundance plot. This is what we see for the microlensed bulge 
dwarfs at [Fe/H] x; -0.4. This means that the star formation 
in the bulge continued for at least as long as it takes for the 
SNIa to contribute significantly to the build-up of elements. 
Estimates of the SNIa time scale depends on the type of envi- 
ronment and rang es from a few hundred mi llion years to 1-2 
billion years ('e.g.. lMatteucci & Recchilf200lh . In Fig. [T] we see 
that the bulge dwarfs with sub-solar [Fe/H] predominantly have 
ages around 10 Gyr, and the lack of an age-metallicity relation 
for the metal-poor dwarfs (Sect. I4.2l i could indicate that SNIa 
started to contribute to the chemical enrichment soon after star 
formation started in the bulge. Even though the internal mea- 
surement errors in the old stars of the sample is still quite large 
to be constraining timescales of 100 Myr versus 1 Gyr, this could 
mean that the timescale for SNIa in the bulge was short. If not, 
the stars with [Fe/H] just above the knee would have been sig- 
nificantly younger than those on the plateau, but they are not. A 
short time-scale of a few hundred million years for SNIa is found 
for elliptical galaxies while the tim e-scale in the disks of spira l 
galaxies are a few billion years ("e.g.. lMatteucci & Recchil200lh . 
Given the similarities we see between the metal-poor bulge and 
the thick disk it is worthwhile to mention that the thick disk in 
the s olar neighbourhood s hows a possible age-metallicity rela - 
tion (IBensbv et al.ll2004al: lHavwoodll2006l: IBensbv et al.ll2007l) . 
with the stars that are more metal-rich than the knee having ages 
that are 1 to 2 Gyr younger than those more metal-poor than the 
position of the knee. Now, we have only 1-2 microlensed dwarf 
stars with -0.4 < [Fe/H] < 0, and it might well be, once we 
have a larger sample, that those stars turn out to be, on average, 
younger than the more metal-poor ones. Also note that there is 
a possibility that the lack of an age-metallicity relation could be 
due to a lot of mixing (either caused by the bar, or caused by the 
merging of sub-clumps). Such mechanisms have been proposed 
to be responsible for both the observed age-velocity relation and 
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Fig. 10. Filled bigger circles show the abundance res ults for 25 
stellar parameters were determined for the star from lEpstein et 
thin disk stars from (jBensbv et al. 2003, 2005i, and in prep.) are 
bars re present the total uncertainty in the abundance ratios and 
dMH, seeSect.E 



of the 26 microlensed dwarf and subgiant stars in the Bulge (only 
ai1 l2010l that we analysed in iBensbv erani20 1 Od) . Thick disk and 



shown as small circles red and blue circles, respectiv ely. The error 
have been calculated according to the prescription in lEpstein et al.l 
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the absenc e of age-metallicity relatio n in the solar neighbour- 
hood (e.g. jMinchev & FamaewOlOl) . 



4.5.3. The metal-rich bulge 

Regarding the metal-rich bulge, most stars seem to follow the 
trends as outlined by the thin disk as seen in the solar neighbour- 
hood. However, there are two stars that have elevated abundance 
ratios that stand out: MOA-2010-BLG-523S at [Fe/H] = +0.09 
andMOA-2010-BLG-259Sat[Fe/H] = +0.37. From Table|2]we 
see that these are two of the more evolved stars in the sample. In 
order to investigate whether the unusually high Na abundances 
in MOA-2010-BLG-523S and MOA-2009-BLG-259S for their 
metallicities could be due to the assumption of LTE, we have 
performed ID non-LTE calculations following lLind et al] ( 1201 lb 
for the four employed Na i lines. We find that while the non-LTE 
effects are significant (-0. 1 dex) they are essentially the same in 
all our targets as well as in the Sun, leaving our derived LTE- 
based [Na/Fe] ratios unaffected. Non-LTE can thus not explain 
the atypical Na abundances of these two stars. Another possi- 
bility could be that the stellar parameters for these stars are er- 
roneous. In order to bring down the [Na/Fe] ratio from +0.19 
to zero for MOA-2010-BLG-523S, and to bring down [Na/Fe] 
ratio from +0.4 to +0.2 for MOA-2009-BLG-259S, we need to 
increase the effective temperatures with more than 500 K or in- 
crease the surface gravities by ~ 1 dex. As the current temper- 
atures of 5250 K and 4953 K well reproduce the wing profiles 
of their Ha lines, it is unlikely that erroneous effective temper- 
atures are responsible for the high abundance ratios that we see 
for some elements in this star Also, both of these stars have more 
than 10 Fe II lines measured and the Fe i-Fe ii ionisation balances 
are well constrained, making a shift of 1 dex in surface gravity 
unlikely. For now we have to accept the fact that these two stars 
show elevated levels for some elements that does not agree with 
the other microlensed bulge dwarf stars at similar [Fe/H]. 

Notable is that the uprising trends of [X/Fe] with [Fe/H] at 
very high metallicities that are seen for the disk stars for many 
elements, e.g., [Na/Fe], [Al/Fe], and in particular the iron-peak 
element [Ni/Fe], is also present in the microlensed bulge dwarfs 
(see Fig. [Toll. 

We furthermore note that some of the characteristics of 
MOA-2009-BLG-259S and MOA-2010-BLG-523S (high Na 
and Al) see m typical of second generation globular clus- 
ter stars (e.g.. ICarretta et al] l2009h . Stars resembling globular 
cluster second generati on stars have been found in the halo 
dMartell & Grebell2010l) and their eventual presence in the bulge 
might gauge the role that globular clusters had, through los- 
ing stars and/or disintegrating in the formation of the bulge. 
However, iMartell & Grebell ( 1201 Oh have not analysed O, Na, 
Mg, and Al, so it is still to be proved whether their sample in- 
deed originated from globular clusters. Also, a connection with 
the globular cluster O-Na anti-correlation for MOA-2009-BLG- 
259S and MOA-2010-BLG-523S is unlikely given that their O, 
Na, Al, and Mg abundances are all high (and several other ele- 
ments too for one of the two stars). 



4.6. Positions and radial velocities 

Figure [TTh shows the positions on the sky for the full sample 
of microlensed dwarf stars in the bulge. They are all located 
between galactic latitudes -2° to -5°, similar to Baade's win- 
dow at {I, b) - (1°, -4°). The reason why most stars are located 
at positive galactic longitudes and negative latitudes is related 
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Fig. 11. Positions, radial velocities, and metallicities for the now 
in total 26 microlensed dwarf and subgiant stars. Grey circles 
mark stars that have [Fe/H] < and black circles those that 
have [Fe/H] > 0. The curved line in (a) shows the outline of the 
southern Bulge base d on observations with the COBE satellite 
dWeiland et alJll994l) . and the concentric dotted lines show the 
angular distance from the Galactic centre in 1° intervals. 



to the extinction / geometry of the bulge. As shown in, e.g., 
iKerins et al.l ( |2009|) , the optical depth, based on stars, is higher 
in these regions and therefore the chances for microlensing to 
occur is higher here. Also, MOA and OGLE have so far mainly 
surveyed the bulge at negative latitudes as these are regions of 
lower extinction. 

We find that there is essentially no difference in the veloc- 
ity distributions between the metal-poor and metal-rich sam- 
ples (see Fig. [TTb). The average velocities are almost identical, 
21.9 ± 98.7 and 19.7 + 109.0 kms"' for the two samples. The 
high velocity dispersion is typical for what has been seen in large 
surveys of the bulge, for instance the BRAVA survey, which ob- 
serv ed and determined ra dial velocities for several thousand gi- 
ants dHoward et al.l2008h . They found velocity dispersions rang- 
ing between 80 - 120 kms"' at these galactic longitudes. 

4. 7. Lithium abundances in the bulge 



As we reported in lBensbv et an(l2010bh . MOA-2010-BLG-285S 
is the first metal-poor dwarf star for which Li has been clearly 
detected in the Galactic bulge. It has a Li abundance which is 
fully consistent with what is seen in other metal-poor dwarf 
stars in the Galactic disk and halo at this effective temperature 
and metallicity (e.g.,|Melendez et al. 2010). Figure [T2h shows 
that the star lies on the metal-rich end of the Li Spite plateau 
(ISpite&Spitelll982l) . Combined with its old age, MOA-2010- 
BLG-285S is an excellent confirmation that the bulge did not 
undergo a large amount of Li production or astration in its earli- 
est phases. 

The Li abundances in the other five more metal-rich dwarf 
stars, in which it could be measured, show a wide range. 
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Fig. 12. (a) Li abundance versus [Fe/H], and (b) versus Teff. 
The large circles indicate the microlensed bulge dwarfs listed 
in Table|6] The ages of the stars have been coded by the sizes of 
the circles a s indicated in (a) . Comp a rison data (small blue dots ) 
con ie from iMelendez et al.1 ( 1201 01) . iLambert & Reddvl (|2004 . 
and lGhezzi et all (l2010albr 



most of them significantly lower than the Li Spite plateau. 
MOA-2008-BLG-311S, which has a higher temperature (^eff = 
5944 K) than the other metal-rich stars, is the one with the 
highest Li abundance. Most likely, the Li in these metal-rich 
stars represent abundances that were obtained by depletion dur- 
ing their main-sequence lifetimes from higher initial values. 
Observations of intermediat e-age solar-metall icity stars both in 
clust ers such as M67 (e.g ., I Jones et all 1 19991) and in the field 
(e.g.. lBaumann et aDl2010h show that they have a large range of 
Li values, which can be explained by a range of ages, masses, 
and initial rotation values. 

The only other detections of Li in the Gal actic bulge are 
from observations of RGB and AGB stars (e.g.. lGonzalez et al] 
l2009l) in which the atmospheric Li abundance have been altered. 
Dwarf and subgiant stars with effective temperatures greater 
than about 5900 K are the only reliable tracers of Li (e.g., 
IWevmann & Searslll965HBoesgaard & Tripiccolll986h . 



does the giant MDF not show this bimodality? Below we discuss 
several possibilities in detail. 

5. 1 . The microlensed dwarf stars are not in the bulge ? 



In iBensbv et al.l ( 1201 Od) we explored in detail whether the mi- 
crolensed dwarf and subgiant stars are located in the Bulge re- 
gion, or in the disk on either this side or the far side of the Bulge. 
We concluded then that the combination of kinematics, colour- 
magnitude diagrams and microlensing statistics makes it highly 
likely that the microlensed dwarfs that we are studying are mem- 
bers of a stellar population that belong in the Bulge. In this study 
we have further investigated the bulge membership of the mi- 
crolensed bulge dwarfs. In Sect. 14. 31 we determined the absolute 
magnitudes of the microlensed dwarfs and compared them to 
those that can be determined from microlensing techniques. We 
find that the two are in agreement at better than our Icr measure- 
ment error of 0.15 mag. This is consistent with the expectation 
(based on Galactic microlensing models) that the overwhelming 
majority of microlensed sources towards these lines of sight are 
in the bulge. 

5.2. Baade's window not representative? 

Could it be that the region towards Baade's window is not rep- 
resentative for the Bulge, not even other close-by regions, and 
that this is the reason for the discrep ant MDFs from the RGB 
and microlensed dwarf star samples? IZoccali et al.l (l2008h find 
different MDFs in their fields (b = -4°, -6°, and -12°) and 
clai m that there i s a ver tical metallicity gradient in the bulge. 
The lZoccali et al.l (J2008') RGB comparison sample is located in 
Baade's window, which indeed is a very small region of the 
Galactic bulge. The 26 microlensed dwarf stars have a wider 
spread in / and b. However, their average distance from the 
Galactic centr e coincides very we ll with Baade's window (see 
Fig. [IB- Also lBrownet al.1 (1201 Ol) used WFC3 photometry and 
two-colour plots to derive MDFs for 4 fields in the bulge in- 
cluding some close to the Galactic plane (b - -2.15°, -2.65°, 
-3.81°, and -4.72°). They found that K-S tests ruled out the pos- 
sibility that they were drawn from the same distribution. They 
were working in reddening-free measurements, and while not 
as good as spectroscopy, obviously, it is a possible confirmation 
of MDF variations across the bulge. Interesting, some of their 
MDFs look quite metal-rich (though they had issues concerning 
their zero-point) and more double-peaked. 



5.3. Strong gradients in the bulge? 

The RGB compa rison sample in Baade's window from 
IZoccali et al.l (l2008h is located 4° from the Galactic plane as well 
as from the Galactic centre, while our microlensed dwarf stars 
range 2°-5° from the Galactic plane and 2°-6° from the Galactic 
centre (see Fig.fTTTi. Would it be possible that there are gradients 
present in the metal-rich bulge population that are not picked up 
by the RGB sample in Baade's window? Possibilities include: 



5. Why are the dwarf and giant MDFs different? 

The MDF is a ke y observable used to constrain chemical evolu - 
tion models (e.g.. lChiappini et alJ[T997t iKobavashi et al.ll2006h . 
Therefore it is important to discuss in detail the differences in 
the MDFs of dwarf and giant stars in the bulge. If the gap in the 
dwarf MDF is real and if the bulge MDF truly is bimodal, why 



- There are two co-existing bulge populations - one metal-poor 
and one metal-rich. However, there is a very steep spatial 
abundance gradient in what we pick up as the metal-rich 
dwarf population, and at the distance of Baade's window the 
[Fe/H] is much lower than close to the plane or close to the 
Galactic centre. Hence, the metal-rich stars can not be seen 
in the RGB sample from Baade's window. 
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- There is a number density gradient in the metal-rich popu- 
lation, meaning that at the distance of Baade's window the 
number of very metal-rich stars has dropped and is very low. 
Hence, the metal-rich population is missing in the RGB sam- 
ple. 



iBabusiaux et al.l ilOKJti use the IZoccaU et all (l2008h data and 
they find a metal-rich and a metal-poor population. They men- 
tion that the metal-rich group seems to drop out as they move 
away from Baade's window, perhaps supporting the second pos- 
sibility wherein the number density of metal-rich star s drops 
quickly. From the same stellar sample iGonzalez et al.l (1201 ih 
show that these two populations separate also in their [a/Fe] 
ratios. However, the fields they are considering are located at 
b - -4°, -6°, and -12°, is a significantly larger range than that 
covered by the microlensed dwarf stars [b = -2°, to -6°, see 
Fig.fTTI) within which the effect should be less dramatic. We also 
note the claim bv iRicb et al. ( 2007i) , which is based on IR spec- 
troscopy, that there is no metallicity gradient in the inner bulge 
(\b\ < 4°). 

5.4. Uncertainties in tlie analysis? 

The analysis of giants is more difficult than that of dwarf stars. 
This is because blending is severe for giants, caused by their 
lower effective temperatures, which will result in spectra with 
stronger lines from neutral elements and a lot of molecular lines. 
The blending will become even more severe for metal-rich stars, 
which could be the reason fo r spuriously hi gh abundance ra- 



tios in some studies (see, e.g.. lLecureuret al.ll2007). Regarding 



i 



S iroble ms in the analysis of solar-metallicity giants. ISantos et all 
20091) reported that their standard analysis of giant stars in open 
clusters results in a much narrower metallicity distribution (with 
a total spread between different clusters of only about 0. 1 dex) 
than that obtained using dwarf stars, which show that there is 
a cluster-to-cluster spread in metallicity about 2.5 times higher 
than that found using giants. The problem seems to be related to 
blends by atomic and molecular lines, as the use of the line list by 
iHekk er & Melendez (2007), which has been checked for blends, 
results in a total metallicity spread derived from op en cluster gi- 
ants s imilar to that derived from dwarfs (Fig. 1 in lSantos et al.l 
|2009|) . Thus, standard analysis of giant stars may result in a com- 
pression of the MDF for metallicities around solar or higher, per- 
haps erasing the metallicity gap we have found in the MDF of 

the bulge. 

The RGB star spectra of IZoccali et al.l (l2008l) have a reso- 
lution of R ~ 20000, whilst our microlensed dwarf stars are 
observed with R ~ 40000 to 60000. That blending might be 
a problem in metal-rich giants can be seen in Fig. [T3h which 
sho ws how the line-to-l ine dispersion of the Fe abundances for 
the IZoccali et al.l ( l2008h RGB sample in Baade's window vary 
with [Fe/H], reaching dispersions of 0.5 dex, and in some cases 
even higher. 

To better account for the large dispersions seen in the RGB 
sample we have in Figure [T3b constructed an alternative MDF 
for the RGB sample. Instead of plotting a standard histogram, 
each star is represented by a normalised Gaussian, defined by 
the specific mean value and error in [Fe/H] for each star. It would 
have been be good if we had error estimates similar to the ones 
that were calculated in Sect.[3]fo r the microlensed dw arfs, but the 
only information we have from lZoccali et alJ (|2P08^ is the line- 
to-line dispersion for each of the 204 stars in Baade's window 
and that in average 60-70 Fe i lines have been used to determine 
the stellar parameters and [Fe/H]. However, Bensby et al. (in 



prep.) have performed an error analysis following JEpstein et al.l 
(1201 Ol) for 703 F and G stars (dwarfs, subgiants, and a hand- 
ful of low-luminosity giants.). It is found that the calculated er- 
rors in [Fe/H] correlate with the line-to-line dispersion of the 
Fei abundances. As the analysis of a majority of the 703 stars 
have u sed many more Fe i lines than the 60-70 that lZoccali et al.l 
(l2008l) have used, we select a subsample of ~ 100 stars from 
the 703 stars that have 100 or less Fei lines measured and used 
in the analysis. We then perform a simple linear regression be- 
tween the line-to-line dispersion and the error in [Fe/H] for those 
stars. This regression line has a constant of 0.044 and a slope of 
0.50, allowing us to transform the line-to-line dispersions into 
estimated errors in [Fe/H] for the 204 RGB stars. This plot is 
sh own in Fig. [T3b. It sh ould be cautioned that the spectra used 
by IZoccali et alj ( l2008l) have lower resolutions and lower S /N 
than the spectra used by Bensby et al. (in prep.), and that the 
stellar parameter determination methodology is di fferent (grav- 
ity is not a spectroscopically tuned param eter in Zoccali et al.l 
l2008l) . Additionaly, and most important, the lZoccaU et alJ(l2008l) 
sample consists of giants while the Bensby et al. (in prep.) have 
analysed dwarf and subgiant stars. However, it is likely that the 
stellar parameters are less well determined when the line-to-line 
scatter is higher. 

The total MDF curve shown in Fig. IT3b is then the sum of 
the Gaussians for all 204 RGB stars, normalised to unit area. It 
is interesting to see how the large uncertainties smears the MDF, 
and especially, now showing a less steep drop-off at the metal- 
rich end. The question now is what happens if the same type 
of uncertainties are added to the sample of microlensed dwarf 
stars? Can it amount for an erasure of the gap such that the dwarf 
MDF resembles the RGB MDF? 

Assigning each star in the microlensed dwarf sample an un- 
certainty given by the dashed line in Fig.lT3b. we then construct a 
normalised gaussian for each of the 26 microlensed dwarf stars. 
These are merged to form a new MDF for the microlensed dwarf 
sample, and is shown in Fig.lT3tl. Although the double-peak still 
is present, it is clear that it has been smeared out. However, the 
double-peak remains. Hence, we can not blame uncertainties in 
the giant analysis to be the sole source for the discrepant MDFs, 
unless the error consists in a compression of the [Fe/H] values 
(see discussion earlier). 

5.5. Disk contamination in RGB sample? 

A certain degree of contamination in the RGB sample in Baade's 
window is expecte d . Based on the Besan con model of the Galaxy 
dRobin et alJl2003h . lZoccaU et alj (l2008l) estimated a contamina- 
tion of 10 - 15 % in their RGB sample in Baade's window, di- 
vided more or less equally between foreground thin disk stars, 
located 2 - 5kpc from the Sun, and thick disk stars located 
within the bulge. The estimation of the thick disk contamination 
is based on the assumption in the Besangon model that the thick 
disk follows an exponential rad ial distribution that p eaks in the 
Galactic centre. As also stated in lZoccali et al.l (l2008l) . these con- 
taminations are poorly determined as we know very little about 
the stellar populations in the inner Galactic disk, and essentially 
nothing about possible thin and thick disk stellar populations 
within the bulge. We therefore find it worthwhile to investigate 
how much contamination is actually needed for the microlensed 

dwarf MDF and the Baade's window RGB M DF to agree. 

T he bulge has a radius of about 1 .5 kpc (iRattenburv et al.l 
l2007l) . meaning that the foreground disk stars, if they are located 
in front of the Bulge, should be located, at most, approximately 
6 kpc from the Sun. Baade's window is located atb - -4°, which 
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Fig. 13. (a) Line-to-line dispersio n of Fe i abundances for the 204 RGB stars from lZoccali et al.l ( l2008h : (b) The estimated random 
errors of the 204 RGB stars from IZoccali et al.l (l2008h . Dashed red line represents the running median (extrapolated outside the 
data limits); (c) The MDF for the 204 RGB stars. The blue line is a smoothed version where each star is represented by a Gaussian 
having with a width equal to the estimated random error as shown in (b); (d) The MDF for the microlensed dwarf stars. The blue line 
is a smoothed version where each star is represented by a Gaussian with a width as given by the dashed line in (b); (e) shows how 
the sum of the minimum residual between a joint disk GCS and microlensed dwarf MDF and the RGB MDF vary with the shift of 
the GCS disk sample. For each value, the fraction of the GCS disk MDF is shown. The minimum sum of residuals is reached when 
the disk GCS MDF is shifted +0.075 dex, and when the fraction of GCS disk MDF is 57 %. The MDFs and the smeared MDFs for 
the RGB sample in Baade's window (f) The smoothed MDFs for the RGB sample (red dash-dotted line), the microlensed dwarf 
sample (green long-dashed line), the GCS201 1 disk sample (olive coloured short-dashed line). The solid blue line shows the mix of 
the MDFs from the GCS201 1 disk sample and the microlensed dwarf sample that best match the RGB MDF. All MDFs have been 
normalised to unit area (hence A^ on the ordinates are arbitrary numbers). 



at this distance from the Sun translates into a vertical distance of 
~ 400 pc. Assuming that the inn er disk contains both a thin and 
a thick disk (see first results in iBensbv et aLluOlOal) . and that 
they have similar scale heights as in the solar neighbourhood, 
this is a distance from the plane where the thin disk is dominat- 
ing. Hence, it is likely that if the giant sample is contaminated 
with foreground stars, it should presumably be by thin disk stars. 

Assigning each star an uncertainty in [Fe/H] given by the 
dashed curve in Fig.[T3b. we create two normalised MDFs: one 
for the bulge (green dashed curve in Fig. [T3f.) and one for the 
disk (olive coloured dashed curve in Fig. [T3f.). The normalised 
bulge MDF is based on the 26 microlensed dwarf stars, and the 
normalised disk MDF on the disk stars from the re-calibration 
of the Geneva-Copenhagen Survey by (ICasagrande et al.ll201 ll 
hereafter GCS2011). In the sol ar neighbourhood the di sk MDF 
peaks slightly below solar (^g.- ICasagrande et al.ll201l1) . As the 
Galactic disk is likely to exhibit a shallow metalli city gradi- 
ent (e.g.. iMaciel & CostallloTol: ICescutti et al.ll2007l and refer- 
ences therein) the disk closer to the bulge would therefore have 
a slightly higher average metallicity. We will therefore leave the 
shift of the GCS201 1 MDF as a free parameter. 

We then construct a joint MDF from the smeared out dwarf 
star MDF and the smeared out GCS201 1 disk MDF, with vary- 
ing fractions of the two (from to 1 in steps of 0.1). Then we 
compare this mixed MDF with the Baade window RGB MDF 
(as given in Fig.[T3h) and determine for which mixture the sum 



of the residuals is minimised. This procedure is done for sev- 
eral different metallicity shifts of the GCS2011 disk MDF The 
results are shown in Fig. [T3fe which shows how the minimum 
sum of the residuals changes with the degree of contamination 
for eleven different shifts of the GCS201 1 MDF. The lowest val- 
ues are achieved for a shift of approximately +0.075 dex of the 
GCS201 1 MDF and a mix of 57 % GCS201 1 and 43 % dwarf 
MDF. As can be seen in Fig. \T3i . the mixed bulge-disk MDF 
and the RGB MDF are similar and agree quite well. Note that 
the exact shift of the GCS201 1 disk MDF seems to be of litde 
importance. All shifts suggest a disk contamination fraction of 
50 to 60 %. 

Is such a high foreground contamination of the RGB sample 
realistic? If 50 to 60% of the bulge RGB samples are actually 
foreground disk stars, this might be visible in the kinematics of 
the RGB sample. However, the velocity dispersion (cr„^) is sim- 
ilar in the RGB sample and in the microlensed dwarf sample. 
Also stars with disk kinematics would introduce a metallicity- 
dependent bias in the RGB sample, since the foreground con- 
tamination only is significant near so lar metallicity (between 
[Fe/H] of -0.3 and +0.1). However. iBabusiaux et al.l (l2OI0h 
found no sign of disk kinematics within that metallicity range. 
Furthermore, the BRAVA survey RGB stars should also presum- 
ably be contaminated by foreground disk stars. However, there 
is good agreement between the iii distribution of the microlensed 
dwarfs and that of the BRAVA giants. 
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Table 7. Properties for the metal-rich and metal-poor bulge sub- 
samples ' . 
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Fig. 14. [Fe/H] versus A ^ax- Grey circ l es mark the new stars and 
black circles the stars in lBensbv et al.l(l2010d) . 



5.6. The magnification puzzle 



In ICohen et al.l (1201 Ol) a peculiar relation between maximum 
magnification (Amax) and [Fe/H] was discovered. Using 16 mi- 
crolensing events available at that time it was demonstrated 
that a very strong correlation between metallicity and maximum 
magnification existed. Several possible explanations were given, 
but all had to be rejected, leaving the puzzle unsolved. Now, with 
additional data points in hand, we again look at this issue. As can 
be seen in Fig. [14] it is evident that most of the events with the 
highest magnifications have super-solar metallicities, and that 
lower magnification events dominate at sub-solar metallicities. 
However, the new data points (grey circles) appears to weaken 
the relation as we now have high-magnification events at low 
metallicities and low-magnification eve nts at high m e tallici ties, 
combinations that were not present in ICohen et alJ (l2010l) . A 
Spearman rank test between [Fe/H] and A^ax now gives a cor- 
relation of 0.42 and a significance that it deviates from zero of 
0.033. A Spearman correlation of zero indicates that there is no 
tendency for [Fe/H] to either increase or decrease when A^ax in- 
creases. This means that we can not rule out the null hypothesis 
that the two properties are uncorrected. So, even though there is 
still a statistical significa nce that the two par ameters are corre- 
lated, it is weaker than in ICohen et al.l (1201 Ol) who found a sig- 
nificance of 0.005, i.e. strongly correlated. More microlensing 
events are needed in order to resolve whether the A„,ax puzzle is 
real or not. 

In any case, assuming that it is real - how would that af- 
fect the microlensed dwarf star MDF? Since all stars observed 
so far, without exception, has turned out to be either metal-poor 
or metal-rich, the only eff'ect would probably be that the rela- 
tive sizes of the two peaks change. The most likely scenario is 
that we are missing more low-magnification events than high- 
magnification events. If so, that would mean that the metal- 
poor peak is under-estimated. The gap around solar metallici- 
ties would probably persist. To investigate what the effects might 
be on the results regarding a possible foreground contamination 
of the RGB sample, we repeat the analysis shown in Fig. [T3f 
with two times as many stars in the metal-poor peak of the mi- 
crolensed dwarf stars MDF. This results in that the best fit is now 
found for a mix of 65 % disk stars from GCS201 1 (now shifted 
0.125 dex) and 35 % microlensed dwarf stars. Thus, if the A,nax 
- [Fe/H] puzzle is real, and should be accounted for, the possi- 
ble foreground contamination of the RGB sample could be even 
higher 





Metal-poor bulge 


Metal-rich bulge 




[Fe/H] < 


[Fe/H] > 
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12 


14 


[Fe/H] 


-0.60 ± 0.30 


0.32 ±0.15 


Age 


11.7±1.7Gyr* 


7.6±3.9Gyr 


[),- 


22 ±99 km S-' 


20 ± 109 kms"' 


[a/Fe] 


0.29 + 0.06 


0.08 ± 0.05 



Notes. ' The values given in the table are the mean values and their associated 
dispersions. Also, note that the ff-element abundance is here defined as the 
mean of the Mg, Si, and Ti abundances. 

+ Excluding MOA-2()10-BLG-446S. If this star is included the average age 
change to 1 1 .2 ± 2.9 Gyr 



5.7. Conclusion 

In summary, we find no definitive answer to why the bulge MDF 
look different depending on if it is based on RGB stars or mi- 
crolensed dwarf stars. Abundance analysis of dwarf stars is cer- 
tainly coupled to smaller errors and uncertainties than that of gi- 
ant stars. Also, the resolution and S /N of the microlensed dwarf 
star spectra are higher than for the RGB stars in Baade's window. 
Furthermore, statistics show that the microlensed dwarf stars are 
highly likely to be located in the Bulge region. However, our 
estimate of the foreground disk contamination is at least four 
times as high as what the Besangon model used bv lZoccali et al.l 
(l2008l) predicts, and is hard to digest. Given that the knowledge 
of the inner disk is very poor, and also the fact that the redden- 
ing is uncertain and often very patchy, it is clear that it is very 
difficult to get a handle on the disk contamination in the RGB 
sample. 

6. The origin of the bulge 

From our studies of 26 microlensed dwarf and subgiant stars we 
find that the bulge MDF is bimodal, and TableQsummarises the 
properties of the two sub-populations. The metal-poor sample 
has an average metallicity of [Fe/H] x -0.60 + 0.30, which is 
in perfect agreement with the mean metallicity th at is found for 
thick disk stars at the solar circle ([Fe/H] - -0.6.^Caro llo et al] 
i201Q) . The metal-rich sample shows a very high-metallicity, and 
narrow, MDF peaking at [Fe/H] = +0.32 + 0.16. The metallic- 
ity dispersion in the metal-poor sample is twice the dispersion in 
the metal-rich one. For the metal-poor bulge we find that almost 
all stars are old with ages in the range 9-13 Gyr, while the stars 
of the metal-rich population have ages that span from the oldest 
stars of the Galactic halo (> 12 Gyr) to the youngest stars of the 
Galactic disk. These age and metallicity differences indicate that 
the bulge consists of two distinct stellar populations; the metal- 
poor bulge and the metal-rich bulge, respectively. Furthermore 
we find that the abundance trends of the metal-poor bulge is 
similar to those of the Galactic thick disk as seen in the solar 
neighbourhood, i.e., enriched [a/Fe] ratios signalling a forma- 
tion on short time-scales. The presence of a "knee" in the [a/Fe] 
abundance plot, coupled with the lack of an age-metallicity re- 
lation, could favour a short time-scale for SNIa in the bulge. 
The similarity between the metal-poor bulge and nearby thick 
disk abundance trends are also seen in s tudies of gia nt stars 
(e.g. Melendez et al. 2008; Al ves-Brito et al . 2010; Rvde et al.l 
1201 Ol) and now this similarity has also been seen between bulge 
giants and inner Galactic disk giants, 3 to 5 kpc from the Sun 
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(iBensbv et al.ll2010a[l201 ih . The abundance trends of the metal- 
rich population partly resemble those of the thin disk in the solar 
neighbourhood. We see no difference in the radial velocity distri- 
butions of the two populations, both show the sam e high velocity 
dispe rsions that is expected for a bulge population jHoward et al.l 
l2009h . 



The observational constraints above points to that the Milky 
Way has a bulge with a complex composition of stellar pop- 
ulations. The similarities between the metal-poor bulge and 
the thick disk suggest that they might be tightly connected. 
Possible scenarios could be that the metal-poor bulge formed 
from an old Galactic disk (secular evolution, evidence for secu- 
lar origins of bulges have been seen in many external galaxies, 
e.g., Genzel et al. 2008), that the metal-poor bulge and the old 
disk formed simultaneously (but separately), or that the metal- 
poor bulge and the thick disk are indeed the same population. 
Simulations have also shown that it is possible for thick disks 
and bulges to form together by rapid internal evolution in un- 
stable, gas-rich, and clumpy disks ( Bournaud et al. 2009) , which 
further could explain the similarities between the Galactic bulge 
and the thick disk. Furthermore, from a differenti al analysis of 
K gian ts at Galactocentric radii 4, 8, and 12 kpc. iBensbv et alj 
(1201 lb found that it is likely that the scale-length of the thick disk 
is shorter than that of the thin disk. This has consequences for 
the relative fractions of thin and thick disk stars in the Galactic 
plane, implying that in the bulge region the two disks should 
be equally present at low galactic latitudes. Interestingly, the 
26 microlensed bulge stars divide about equally into being ei- 
ther metal-poor or metal-rich (see Fig.|6l). Could this mean that 
what we are seeing is a bulge solely made up of thin and thick 
disk stars, and that the metal-rich peak in the microlensed bulge 
dwarf MDF is a manifestation of the inner thin disk? The thin 
disk is known to sh ow a strong metallicity gradient (see, e.g., 
ICescutti et al.l 120071) . and if interpolated (linearly) to the bulge 
region it would reach [Fe/H] > -1-0.5. A linear relationship all the 
way to the bulge is unlikely, but the strong gradient indicates that 
the average metallicity could end up at, or around, the metalli- 
icty that we see for the metal-rich peak of the microlensed dwarf 
MDF ([Fe/H] » -1-0.3). The lack of a radial gradient in the thick 
disk could be explained through radial migration that has had 
time to wash out abundance gradients in this old population. For 
the much younger thin disk, radial migration has not had enough 
time to act, and a radial metallicity gradient is still present. These 
findings could support the results from the BRAVA su rvey which 
claim that there is no evidence of a classical bulge dShen et al.l 
I2OIOI) . but that the Milky Way is a pure-disk galaxy. The forma- 
tion of the bulge f rom a two-component stell ar disk has recently 



been modelled b y Bekki & Tsujimota (|201 ll) . We note that also 
iBabusiaux et al.l (1201 Oh find evidence for two bulg e populations. 
Their results are based on the metallicities from IZoccali et al.l 
(I2OO8I) . now complemented with kinematics, and they find that 
the metallicity grad ient along the bulge minor axis observed by 
IZoccali et al.l (l2008h can be related to a varying mix a metal-rich 
population with average metallicity of [Fe/H] = +0.13 that has 
bar-like kinematics and a metal-poor population with average 
metallicity of [Fe/H] - -0.3 that has thick disk-like kinematics. 
Regarding the metal-rich bulge we find a wide ra nge of stellar 
ages. As also discussed in IBabusiaux et al.l (1201 Ol) . this is con- 
sistent with a stellar population with old stars from the inner 
disk that have been re-distributed by the central bar, and whose 
young stars have formed in star bursts triggered by the central 
bar. 



7. Implications for the bulge IMF 

The MDFs based on the microlensed bulge dwarfs and on red 
giants in Baades window are clearly different. This has impor- 
tant consequences for chemical evolution models, as the MDF 
strictly constrains which kind of models are plausible. In par- 
ticular, the peak of the MDF depends on the slope of the initial 
m ass function (IMF). Th e chemical evolution model of the bulge 
bv lBaller o et al.' ('2007b') is based on the photometric MDF of gi- 
ant stars by Zocc ali et al. (2003) and on t he spectroscop ic MDF 
of giant stars by iFulbright et al] ( |2006[) . iBallero et all ( l2007bl) 
suggested for the bulge an IMF much flatter than in the solar 
neighbourhood, i.e., an IMF skewed towards high mass s tars, 
to explai n the MDF o f bulge giants by IZoccali et alj (l2003h and 
Fulbright et al.l (|2006|) . 

A more recent model bv ICescutti & Matteucci (1201 ll). us- 
ing no w the spectroscopic MDF of giant stars by IZoccali et al.l 
( I2OO8I) . also favours a flat IMF to reproduce the peak of the MDF 
of bulge giants. However, those conclusions may have to be re- 
vised in view of our new MDF based on microlensed dwarf stars. 
Our MDF has two well-defined peaks, one that we associate with 
a metal-poor old bulge, and another one at super-solar metallic- 
ities associated with a younger population. Hence, the IMF no 
longer has to be flat to explain a single p eaked solar metallicity 
MDF that was made in 0.5 Gyr, as in the ICescutti & Matteuccil 
(l2011h models. To reproduce the bimodal nature of the bulge 
MDF, probably a normal IMF can be used for the metal-poor 
bulge, while contributions from Type la SN can probably ex- 
plain the younger metal-rich peak. 

Notice that our MDF is still poorly sampled (only about two 
dwarfs per metallicity bin), thus, in order to provide even firmer 
constraints on chemical evolution models, it is important to at 
least double the sample of microlensed bulge stars observed at 
high resolution and high S IN. This will provide tight constraints 
on our galactic bulge, and on a broader context will also have 
impo rtant consequences o n the modelling of external galaxies 
(e.g.. lBalleroetal.ll2007ah 

8. Summary 

With the twelve events presented in this study the sample size of 
microlensed dwarf and subgiant stars in the Galactic bulge has 
grown to 26. All stars have been observed with high-resolution 
spectrographs on 8-10 m class telescopes, which has allowed us 
to determine detailed elemental abundances for many elements. 
The following results and conclusions can be drawn with the 
current sample. 

1 . The MDF based on the microlensed dwarf and subgiant stars 
is bimodal with one peak at [Fe/H] ^ -0.6, one peak at 
[Fe/H] X! -1-0.3, and with essentially no stars around solar 
[Fe/H] . This is in stark contrast to the MDF provided by giant 
stars in Baade's window which peaks at this exact location. 

2. A statistical test shows that it is unlikely to have a gap in 
the MDF based on the microlensed dwarf st ars if the mi- 
crolen sed sample and the giant sample from IZoccali et al] 
(I2OO8I) have been drawn from the same underlying popu- 
lation. We find no definitive answer to why the two distri- 
butions look different. Possible explanations include: large 
uncertainties in the analysis of the giant stars (especially the 
metal-rich ones); a higher level than predicted contamination 
of foreground disk giant stars in the giant sample (maybe as 
high as 60 %); low number of microlensed dwarf stars, the 
difference will go away once more microlensing events have 
been observed; a mix of all above. 
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3. The metal-poor bulge is very similar to the Galactic thick 
disk in terms of mean metallicity, elemental abundance 
trends, as well as stellar ages. The metal-rich bulge is com- 
plicated, resembling the thin disk in terms of abundance ra- 
tios, but the thin disk does not show a wide range in ages. 
These findings could be evidence of the co-existence of 
two distinct stellar populations within the bulge, each which 
might have its own formation scenario. 

4. Combining the similarities between the metal-poor 
bulge and the thick dis k with add it ional evidence from 
iBensbv et all (1201 lb and IShen et al.l (1201 Oh . we speculate 
that the bulge is purely made from the stellar disk, where the 
metal -rich bulge would be the manifestation of the inner thin 
disk. The origin of the metal-rich bulge population might 
also be coupled to the Galactic bar, including old stars from 
the inner disk scattered by the bar and younger stars whose 
formation was triggered by the bar. 

5. Li abundances were presented for six of the microlensed 
stars. One star is located on the Spite plateau ([Fe/H] - 
- 1 .23 and log e(Li) = 2.16), indicating that the Spite plateau 
might be universal. The other five stars are all metal-rich and 
show evidence of Li depletion due to their lower effective 
temperatures. 

6. The new events presented in this study have weakened the 
recently discovered, and un-explained, Ajnax-[Fe/H] puzzle. 
We now have stars at high [Fe/H] that have low A^ax and 
stars at low [Fe/H] that have low Amax- More events are 
needed to reveal the t rue nature of this puzz l ing re lation. 

7. The recent claims bv lCescutti & Matteuccil (1201 lb of a very 
flat IMF in the bulge is not supported by the MDF and abun- 
dance trends as probed by our microlensed dwarf stars. 

8. Lastly, using the new colour-temperature calibration by 
ICasagrande et alj ( 1201 Ol) and the now in total 26 microlens- 
ing events, we find that the colour of the bulge red clump is 
(y-/)o = L06. 
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